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ABSTRACT 

Gamma ray line emission processes in solar flares and cosmic 
transients are reviewed and implications of recent line observations 
are discussed. The gamma ray line emission from solar flares results 
from nuclear interactions of accelerated particles with the solar 
atmosphere. The observed line intensities give information on the 
total number and spectrum of particles accelerated in the flare , on 
the temperature and density in the interaction regions and on the 
time history of the interactions. Analysis of the line observations 
from the June 7, 1980, flare show that the number of protons acceler- 
ated in the flare exceeded the number observed in the interplanetary 
medium by a factor of -100. The bulk of the accelerated protons, 
therefore, remained trapped in the solar atmosphere. where they pro- 
duced gamma ray line emission as they slowed down. 

The gannr-i ray emission lines observed from gamma ray transieiits 
appear to result from both the annihilation of positrons produced in 
photon-photon interactions and from deexcitation of nuclear levels 
and capture of neutrons produced in nuclear interactions. The observ- 
ed line intensity provides information on the temperature, density 
composition, magnetic field and redshift in the transient sources. 

Both the gravitational reashift and the iron enrichment implied by the 
gamma ray line observations strongly suggest that neutron stars are 
the soxirce of the trsuisients. 

•Research supported by NASA grant NSG-75^1- 
••Also at University of Maryland, College Park, MD., research sup- 
ported in part by NASA’s Solar Terrestrial Theory Program - 
Permanent address: Tel Aviv University ^ Ramat Aviv Israel, 
research supported in part by the U.S. -Israel Binational Science 
Foundation. 
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INTRODUCTION 

Gamma ray lines are unique tracers of the high energy processes 
that dominate the physics of solar flares and cosmic transients. In 
solar flares the lines are produced by nuclear interactions of ener- 
getic protons and nuclei with the solar atmosphere. Similar inter- 
actions in high temperature plasmas may also produce geunma ray lines 
in cosmic transients. In addition, the extremely high photon 
densities, matter densities and magnetic fields expected in the 
sources of these cosr.dc transients should lead to gamma ray line pro- 
duction processes that have no counterparts in solar flares. Because 
of the high photon densities, photon-photon collisions can become a 
dominant line producing mechanism via e^-e” pair production. The 
high pair densities resulting from these interactions may also lead 
to pair degeneracy and to possible maser action. The ultra high 
magnetic fields produce cyclotron emission and absorption lines, end 
cool and confine the e"*‘-e“ plasma, a necessary ingredient for pro- 
ducing an observable smnihilation line. 

The strongest line in solar flares is at 2 .t ^23 MeV from neutron 
capture on hydrogen in the photosphere. This result, predicted by 
theory' and now confirmed by many observations, is consistent with 
line production by accelerated particle interactions. Indeed, for 
particles with essentially cosmic abundances and energies around a 
few tens of MeV/nucleon, neutron productioxi greatly exceeds the pro- 
duction of any single nuclear deexcitation line. The prompt deexci- 
tation lines, such as the MeV line of do, however, provide 

a direct measure of the energetic particle reaction rate. This line 
as well as other prompt lines have been observed from several flares, 
revealing the time history of the nuclear interactions. 

The strongest emission line, observed in gamma ray transients, is 
that in the energy range from O.i^O to O.I46 MeV. It is reasonable to 
associate this line with e'''-e“ annihilation. The redshift from a 
line energy of O.51I MeV or more could be due to the gravitational 
redshift of a neutron star or to maser (grasar) action discussed eise- 
where in this volume. Lines from nuclear reactions also seem to be 
present in cosmic transients. In addition to the gravitational red- 
shift, the principal difference between these lines and the solar 
lines is the drastically different particle abundance needed to 
explain the cosmic transient lines. In particular, an enhanced 
abundance of ^^e is required to account for both a nuclear deexci- 
tation line in a gamma ray burst and a neutron capture feature in a 
longer duration transient. This result also supports a neutron star 
origin for cosmic transients. 

In the first part of this paper we review the theory of solar 
gamma ray line production and present results of new numerical 
calculations based on more detailed and accurate nuclear cross 
sections. We illustrate the application of the theory by considering 
the June T, I98O, flare for which there are reasonably detailed 
gamma ray data from SMM and important supporting particle data from 
ISEE- 3 , Helios -1 and IMP - 8 . 
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In the second part of the paper we discuss the physical processes 
responsible for ganuna ray line production in cosmic tremsients, in 
particular, positron production and annihilation emd nuclear line 
emission. We also review the possible origin of the gamma ray lines 
observed in a longer duration transient which is a particularly strong 
gamma ray line emitter. 


SOLAR FLARES 

Nuclear interactions of flare-accelerated protons and nuclei in 
the solar atmosphere p'coduce gamma ray lines from neutron capture, 
positron annihilation and nuclear deexcitation. The first solar 
, gamma ray lines at 0.511, 2.22, and 6.13 MeV, wore observed 

by Chupp et al^ with a Nal detente' on the OSO-7 satellite during the 
August U, 1972 flare. These an >ther lines nave been observed from 
several subsequent flares by detectors on the HEAO-1^, KEkO-3^ and 
SMM** satellites. The relative intensities of these lines were con- 
sistent with earlier predictions.^ A review of the observations is 
given by Chupp in this volume. 

Neutron, Positron and Deexcitation Line Production 

The strongest gamma ray line observed in nearly all of th^se flares 
is that at 2.223 MeV from neutron capture on hydrogen, '^Ii(n,Y) n. 

Several theoretical studies 5-10 been made of neutron production, 

slowing down and capture in solar flares. Here we briefly review the 
interaction models used and we present results of new numerical 
calculations based on updated nuclear cross sections, which will be 
published elsewhere. 

The ii;stantaneous, or tliin-target , neutron production rate is 
given by 

«BOjj(E)Nj(E) (1) 

where n^. is the ambient hydrogen density, a^ is the relative number 
of target nuclei i with respect to hydrogen, and 0 . . is the neutron 
production cross section from the interaction of tai^get nuclei i with 
accelerated particles J having velocity 6c, energy per nucleon E, and 
energy dependent number N.(E). 

^ Ihe time-integrated, ^or total thick-target, neutron yield, is 

given by 


where A^ is Avogadro's number, (dE/dx)j is the energy loss rate per 
nucleon of accelerated particles J per^^xinit path length of the ambient 
medium, and N*(E) is the total number of accelerated particles J 
incident on tne target. For more details see ref. (lO). 

Various for^ of accelerated particle spectra were used in pre- 
vious treatments^ ~ of neutron and gamma ray line production ^n 
energetic particle reactions. Two such fo:*ris, which also give^ 
good fits to the particle data in interplanetary space, are exponen- 
tials ir. rigidity 
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Nj(E) or iJjCE) = Aj exp(- Pj/P^) dP^/dE (3) 


and Bessel functions, 

Nj(E) or Nj(E) = Aj K 2 ^ 2 ( Sp/m^CQ T)] 


(ll) 


Here A, in proportional to the abundance o f the enerfi etic pirticles 
J, P, = (A/Z) p is particle rigidity, p = /E(F + Pm^c^^yis particle 
momentum per nucleon, m is the proton mass. K 2 is the modified 
Bessel function of orde? ?, and P^(measured in MV) and a T(nondi- 
mensional) characterize the spectrum of the energetic particles. 

Equation {I 4 ) was shown^^ to be the particle number spectrum that 
results from stochastic Fermi acceleration with no energy losses, an 
energy and charge independent acceleration efficiency a and an escape 
time T. Furthermore, for the June 7, 1980, flare which is the first 
SMM event with reasonably complete gamma ray line and interplanetary 
particle data, this Bessel function spectrum gives a very good fit to 
both the proton and alpha particle spectra with essentially the same 


aT. 

The solar atmospheric abundances^^ and the energetic particle 
abundances, used in the present calculations, are listed in Table 1. 
The particle abundances, which show a substantial enrichment of heavy 
nuclei, are known (e.g. ref, I 6 ) to vary from flare to flare. 

The cross sections for neutron production in energetic particle 
reactions have been summarized in ref, (lO), We have recently up- 
dated these cross sections and added many new reactions that involve 
all the isotopes listed in Table 1, These new cross sections, to be 
published elsewhere, are used in the present calculations. 

The magnitude, spectral form and charge dependence of the energy 
loss rate (dE/dx)j have important consequences on the nuclear yields 
in thick targets. 

The present calculations have been carried out for slowing down 
in a neutral medium. The expression 


(dE/dx) = (Z^ ff/A). 630E"°‘® (MeV/nucleon)(g/cm^)”^ (5) 

J ^ J 


where 

= Z [l-exp(-13TB/Z^^^)] , ( 6 ) 

IT 

gives a good fit to the tabulated values of dE/dx for charged par- 
ticles slowing down in neutral H, In an ionized medium the energy 
losses are higher^^ and hence the gamma ray yield per unit energy 
deposited is lower. We limit the present treatment, however, to a 
neutral medium. 
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tabu; 1 Abundances 



Ambient 

Energetic 

Isotope 

Particles 

Particles 

h 

1. 


1. 


**He 

0.07 


0.15 



It. 15 X 

10“'* 

1.07 

X 10"3 


lt.6li X 

10 "^ 

1.28 

X 10"^ 


9.0 X 10“^ 

2. lit 

X lO"** 

'5n 

3.*t6 X 

10"^ 

8.57 

X 10"^ 


6.92 X 

10-*^ 

2. lit 

X 10"3 


1.38 X 

10"^ 

It. 28 

X 10"^ 

2 °Ne 

9.0 X 10"^ 

2. lit 

X 10"** 

^"Ne 

1.0 X 10"^ 

2.57 

X 10"^ 

23„ 

Na 

P.28 X 

10 -^ 

It. 28 

- 10"^ 


3.11 X 

10 “^ 

6.1t2 

X 10"** 


It. 01 X 

10 -^ 

8. lit 

X 10"^ 


It .lt3 X 

10 -^ 

8.it9 

X 10"^ 


3 .i 8 X 

10 -^ 

5.35 

X IC"^ 

"®Si 

3.>t6 X 

10 "^ 

6.1t2 

X 10"** 

?9si 

1.60 ■ 

10 -^ 

3.21 

X 10"^ 

3°Si 

l.l8 X 

10-^ 

2. lit 

X 10"^ 

32s 

1.80 X 

10“^ 

1.07 

X 10"** 

3^ 

7.61 X 

10 

It. 71 

X 10"^ 

3“at 

3.39 X 

10-^ 

2. lit 

X 10"^ 

36Ar 

6.23 X 

10 “^ 

It. 28 

X 10“^ 

“^CS 

2.28 X 

10-^ 

It. 28 

X 10"^ 

52cr 

Ji.l5 X 

10" 

2. lit 

X 10"^ 


I.9I1 X 

10-^ 

6.85 

X 10"^ 


3.11 X 

10-5 

1.07 

X 10-3 

5-^>-e 

7.61 X 

lo"-^ 

2.57 

X 10"^ 


1.25 X 

10-^ 

2. lit 

X 10"^ 


1 ,8tt X 

10"’^ 

8.57 

X 10"^ 
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Fig. 1. Total neutron production in the thick- and thin-target 
interaction models with Bessel function (oT) and exponential 
(Pq) spectra. 

The calculated neutron production rates and yields are sha^ in 
Figure 1 for energetic particle densities normalized to i proton above 
30 MeV. For relatively flat energetic particle spectra, corresponding 
to large values of or oT, the bulk of the neutrons are produced in 
pa reactions. For steep particle spectra, given by the Bessel fiincUon 
at small aT, the large neutron yields result from reactions between 
a-particles and heavy nuclei. This effect is absent for the rigidity 
spectra because particles with Zi2 have lower energies per nucleon at 
the same rigidity and hence produce less nuclear reactions. These 
effects were discussed in more detail in ref. (10) in connection with 
the comparison cf neutron production by particles with spectra that 
were either power laws in energy or exponentials in rigiuity. 

The 2.223 MeV geunma ray line emissivity is equal to the neutron 
production rate times a queintity f which is the product of the 
fraction of neutrons that are captured on hydrogen in the photosphere 
and the probability of the resultant 2.223 MeV photon escaping from 
the photosphere. These neutron captures must compete with nonradiative 
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captures, as well as with neutron decay ai;d escape from the svin. Non- 
radiative capture on ^He is the most important competing reaction®, 
even though ^He is only a minor constituent of the photosj^ere, 
because the cross section for the reaction ®He (n,p)®H is about four 
orders of magnitude larger than that for the reaction (n,y)^H. Ob- 
servations of the intensity of the 2.223 MeV line compared to that of 
other lines can limit the photospheric ®He/**He ratio. Calculations 

asrumlng ^He/^H of 5x10“^, isotropic neutron production above 
the photosphere and a flare site away frcm the limb, give ? ranging 
from 0.1 to O.ll*, depending on the energy spectrum and interaction 
model. In the thick target case, however, the neutrons could be 
produced in the photosphere which could increase ? to as much as ”0.2. 


The 0.511 MeV line from positron annihilation has been observed 
from several flares. A number of theoretical studies have been made 
of positron production 5,7,10 positron slowing down and annihi- 
lation For this paper we give the res\ilts of new calculations 

of positron production based on much greater number of 0+ emitters 
than were considered in the previous calculations. The results are 
given in Figure 2, where we show the ratio for thick and thin targets, 
Qe+/Qn 'le+Z^n respectively. 


oT 



Figure 2. Ratios of the totfj. position yield to the total neutron 
yield for the thick -and thin-target interaction models and Bessel 
function (oT) and exponential (Pq) spectra. 

and were calculated using the same equations, abundances and 
particle spectra as given above for the neutron calculation. The 
positron yields shown in this figure represent total yields. Because 
of the finite half-lives of the various 6+ emitters in a short obser- 
vation time of a transient event, fewer positrons than indicated in 
Figure 2 are available for 0.511 MeV line production. 
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Flare-accelerated particle interactions also lead to many other 
camma ray lines from deexcitation of nuclear levels. The two strong- 
est lines, at and 6.13 MeV due to and ^^0* deexcitation 

respectively, were first observed* from the solar flare of August 
197<^* These and other nuclear l.'nes have beer seen in a number of 
subsequent flares, reviewed by Chupp elsewhere in this volume. 

We have treated*^ in detail the production of gamma ray lines in 
energetic particle reactions. Using this treatment, we have evaluated 
the prompt gamma ray lines fbr the abundances, particle spectra and 
interaction models discussed above. Each of the resultant lines has 
two components: a narrow (aE/E = 2,^f) line component from deexcita- 
tion of the lecoiling, ambient-gas nuclei emd a much broader under- 
lying component from deexcitation of the fast, accelerated-particle 
nuclei. The bxilk of the gamma ray flux, observed from the August *4, 

1972 flare at energies between U and 7 MeV, has been shown** to result • 

from the superposition of broad and narrow nuclear lines rather than 

electron tremsstrahiung or other continuum emission processes. This 

energy band, now referred to in the SMM observations as the "main 

channel window", can thus provide a direct and sensitive measure of 

this interaction rate of flare accelerated nuclei. 

Because the cross sections of nuclear excitation levels have 
different energy dependences from those of neutron and positron pro- 
duction, the calculated ^**® ratios of nuclear deexcitation line 
emissivities to the neutron production rate depend strongly on the 
assumed spectra of the accelerated particles. One of the strongest 
nuclear deexcitation lines is that at MeV. 




«T pj(MV) 


Figures 3 and 1*. Patios of the I 4 .UI 4 MeV to the 2.22 MeV line 
intensity in the thick-and thin-target models with f = 0.12 for 
Bessel function (aT) and exponential (Pq) energetic particle spectra. 
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In Figures 3 and ve show ratios of the intensity of this line 
to the 2.223 MeV line intensity for the different interaction models 
and energetic particle spectra. As can be seen* these ratios generally 
decrease with increasing energetic particle spectral hardness reflect- 
ing an increased neutron production and decreased MeV photon 
production by particles of high energies. Likewise ^(l4.l4l*)/^(2.22) 
is lower for thick targets than for thin targets (except for very 
steep Bessel function spectra) because the energy losses harden the 
particle spectra in the thick target. For very steep (smalZ oT) 

Bessel function spectra, the thick target ratio exceeds that for thin 
targets, because of the suppression of the heavy particle fluxes by 
energy losses in the thick target. As mentioned above, the inter- 
action of these particles with helium are the main neutron sources 
at low energies. This effect is not present for rigidity spectra, 
because, as explained earlier, the contributions of the heavy particles 
to the nuclear reaction rates is negligible for such spectra for both 
thick and thin targets. 

Implications of Line Observations 

The observatioi.d of the intensity, time history and relative 
strength of gamma ray lines from solar flares allow us to determine 
important properties of the flare accelerated particles, the interac- 
tion and emission region and the particle acceleration process. 
Particle Spectrum and Number 

In this paper we shall focus on the flare of June 7, I960 for 
which there are both gamma ray line observations and interplanetary 
particle measurements. The combined analysis of these observations 
and measurements provide important new information. In particular 
we find that in this flare the bulk of the gamma ray line emission 
results from thick-target interactions of accelerated particles 
trapped and slowing down in the solar atmosphere and not from inter- 
actions of the particles during their acceleration. Moreover, for 
this flare the gamma ray line intensities require that the bulk of the 
flare accelerated particles remained trapped in the solar atmosphere 
and that only a very small fraction of the accelerated particles 
escaped into the Interplauietary medium. 

The June 7 flare was located at N12W7U and hence was well 
connected magnetically to Earth. Indeed, several observations of 
energetic particles have been reported^ ^ . Based on these, the 
number of protons of energies greater than 10 MeV released into the 
interplanetary medium, Nesc( >10MeV ) , has been estimated^ ^ Lo be '10^*. 
Furthermoie, the spectrum of these protons was well fit with the 
Bessell function of equation {h) with a? equal to 0.013 (R. McGuire, 
private communication 198 I). This spectral form also fits the a* 
particle spectrum with essentially the same oT. 

1.-. the thin target model, the 2.22 MeV line fluence can be 
written^ ^ as 

(2.22) « 7 MeV)/i*nd2 (7) 

here T is the particle escape time from the thin-target interaction 
region and d is the distance to the sun. If this region is also the 
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acceleration reelon, the time T in equation (7) is the sane as the 
time T in the parameter oT that defines the spectral hardness of the 
Bessell function. 

Using the proton data (aT*0.013, and Nesc(>10MeV) ■10^* protons 
giving Nesc(>3oMeV) ' 5x10^*) and the gamma ray data (^(2.22) “6.6 
photons cTn“^from Chupp in this volume), we obtain from Figure 1 and 
equation (7) that for a thin target n^T “ 7 x 10*** cm^ sec”*, or a 
matter traversal for 30MeV crotons of 12p;nicni“2 . The large abundanca? 
spallation products (^H, Li, Be, B) that would result from such a 
long path length are not observed in solar flares. We conclude, 
therefore, that the gamma ray lines observed from the June 7, I980 
flare were probably not produced in thin-target interactions. A 
si.nilar conclusion has been reached*^ for the August *4, 1972 flare. 

In the thick-target models, on the other hand, the spallation 
products, that accompany the production of gamma ray lines, are slowed 
down in the solar atmosphere and hence are not expected to be seen in 
the interplanetary medium. In this model the 2,22 MeV fluenc . is 
given by 

^(2.22)-?QnRp(>30MeV)/i47Td2 (8) 


where 0^ is calculated from equation (2), 

If we assume that the same process accelerates both the particles, 
!Ip, that interact at the sun to produce the gamma ray lines and the 
particles, Nesc, that are observed in the interplanetary medium, and 
that both populations therefore have essentially the same energy 
spectrum, then Rp should be appioximately equal to Mesc (l-F)/F, 
where F is the fraction of particles escaping into the Interplanetary 
medium. From the calculated for aT'0.013 (Figure 1) and the 
ol'served (^(2.22)* 6.6 photons cm”^ we then find that Rp (>30MeV) * 
14.5x10^* protons or Rp (>10MeV)= 10^^ protons. This exceeds 
Nesc(>10MeV) by a factor of '100, implying that F-10"^ or only 'ij 
of the accelerated protons escaped from the June 7 flare and that the 
rest remained trapped in the solar atmosphere where they produced the 
observed gamma ray line emission by nuclear interactions while they 
slowed down. 


Analysis of the prompt gamma ray line emission from the June T 
flare also suggest that the observed emission was produced in thick- 
target interactions, not thin target. In particular, we consider the 
ratio of the combined line fluerce in the ^4 to T MeV band to the 
fluence in the 2.22 MeV line, which in this flare was found to be 
1.7l4r0.27 (from Chupp in this volume). 

Using the methods and nuclear data oi ref. (12) with the Bessel 


function spectrxim of equation (-4) and the abundances given in Table 
1, we calculate for aT«0.013 that q(i4-7)/nH»!4.i4xlO“**’Np( >30MeV) erc^ 
se'.'* for the thin-target model and Q( l4*-7)«l ,ixl0"^Np( >30MeV) for the 
thick-target model. Comparing with the neutron production rates 
(Figure 1), taking f ' 0.12, we calculate ; '-T) /^(2,22) = 2.6 for the 


thick target and "7.5 for thin target. The observed ratio of 1.7*4* 


0.27 is in much better agreement with the thick-target model than with 


the thin-target. Moreover, since in the thick-target model the 
ne-trons ma>’ be produced in the photosphere where ? could be as much 
as 'C.2, the calculated $( *4-7 ) /♦( 2 . 22 ) could be as low as '1.6 and 


« 
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thus be in excellent acreement with the observations. Other effects, 
such as some beaming of the charged particles, could also produce 
similar agreement in the thick^target case. 

St ructure of the Interactio n Region 

Theoretical studies 5,0-10 qj solar flare gamma ray line 
emission processes suggest that the prompt, deexcitation line 
emission, the subsequent positron annihilation and the evertual 
neutron capture emission should each occur at successively deeper 
mean depths in the solar atmosphere. This results from two effects. 

* First, the range of the 'lOMeV protons primarily responsible for the 
excitation of nuclear levels is only about l/5th of that of the 
*30MeV protons that produce the bulk of the positrons and neutrons. 
Thus, if these particles are accelerated high in the solar atmosphere 

* and travel down the magnetic field lines into the deeper atmosphere 
where they lose their energy and are eventually stopped, the bulk of 
the nuclear line excitation should occur at a higher altitude than 
the neutron and positron production. Second, since the range of the 
secondary positrons is less than that of the "30 MeV protons which 
produced them, the positrons should annihilate at a depth close to 
that at which they were produced. But the neutrons on the other hand 
have a slowing down and capture mean free path that is greater 

than that of their proton pregenitors, so they can propagate to sig- 
nificantly greater depth before they are captured. On the average, 
therefore, those neutrons that are capti.red will do so at a greater 
average depth than that at which they were produced, since those 
neutrons which move upward to shallower depths are more likely to 
escape or decay before they can be captured. 

Strong observational evidence for such differences in the ,aean 
emisrion depths of the deexcitation and capture lines appeari. in the 
gamma ray line observations of the limb flares of June 21, 1500 and 
April 27, 1901, described by Chupp in this volume. In these . .ires 
the observed ratio of the deexcitation line fluence at h-7 MeV to that 
of the 2. 22M'iV neutron capture line wac roughly a factor of times 
and 10 times that of the average disk flares, respectively. 

This is clear observational eviaence that the 2.22 leV line is 
produced in the photosphere with differential attenuation, or limb 
darkening, of the 2.22 MeV line corresponding to column density 
differences between the nuclear deexcitation region and the neutron 

^ capture region of (2 to 3) x 10 ^ cm"^. The exact longitudes of these 

limb flares are not known but these differences probably correspond 
to radial column depth differences of a few times cm“^. S\xch a 

value for the depth of the mean neutron canture region is also 

* consistent with the average den.^ity of cm”^ in the neutron 

capture region for the June 7, I9OO flare implied by the 'xO^ sec 
decay tine of its 2.2 MeV line emission. Detailed studies of the 
time dependence of the positron annihilation and neutron capture line 
emission in all of these flares should provide much more information 
on the structure of the flare emission region. 
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GAI-IMA RAY BURSTS 

Observations of emission lines and absorption features, seen in 
the spectra of many gamma ray bursts and transients, are reviewed 
by Teegarden in this volume. 

The absorption features are observed^** at energies below about 
100 keV in the spectra of 20 bursts. If tnese are due to cyclotron 
absorption, they require very strong magnetic fields of the order of 
10^^ gauss, such as those expected around neutron stars. The most 
commonly observed emission line falls in the range from O.U to O.I 46 
MeV, as observed^** by low resolution Nal detectors in seven gamma ray 
bursts. In the spectrum of one of these bursts, that of November I 9 , 
1978 , a small Ge detector has resolved^ ^ two lines at ~ 0.i*2 MeV and 
■"O.T^ MeV, which the Nal detectors have seen as one broad emission 
feature from 0.3 to 0.8 MeV. Line emission in the range of O.U to 
0.i»6 MeV is most likely due to gravitationally redshifted e'*‘-e“ 
annihilation radiation, while the line at 0.7^ MeV could be either 
collisionally excited and gravitationally redshifted 0.8i*7 MeV 
emission from ^^Fe, an abundant constituent of the crusts of neutron 
stars, or gravitationally redshifted single photon e'*'-e“ annihila- 
tion^^>^^ radiation at 1.022 MeV in a very strong ("lO^^ gauss) 
magnetic field. In all cases, the implied redshifts of 0.1 to 0.3* 
are consistent with those expected from neutron star surfaces. Thus, 
the magnetic fields, the redshifts and the surface composition indi- 
cate that the sources of the bursts with observed lines are probably 
neutron stars. 


Positron Annihilation 


The principal positron producing mechanism in gamma ray b\irsts is 
likely to be pair production in photon-photon collisions. This follows 
from the enormous photon densities in the burst sources deduced^® 
from the observed liuninosities , likely source distances and sizes. The 
latter probably are of the dimensions of neutron stars or smaller. 

The production of a relatively narrow e"‘’-e“ annihilation line 
req'jires®® that the pairs annihilate at a temperat\ire v'hich is sv.b- 
stantially lower than the burst temperature deduced from the observed 
continuum spectrum or the temperature required to produce the pairs. 

A promising cooling mechanism®^ is synchrotron radiation in strong 
magnetic fields ( B>10^ ^gauss ) . The radiation produced by this cooling 
could, be responsible®^*®^ for the continuum emission of the bursts at 
energies below -300 keV. 

This synchrotron cooling model was specifically applied ®^>®2 to 
the March f>, 1979 transient. Although the observation of a '0.l43 MeV 
emission line from this transient suggests that its source was a 
neutron star, the origin of this spectacular event remains unresolved. 
The position of the burst has been determined by triangulation (see 
review by Cline in this volume) leading to an error box of size less 
than an arc minute within the supernova remnant N^*9 in the Large 
Magellanic Cloud (LMC). Nevertheless, the question remains whether 
the burst did indeed originate in the LMC or whether its source was 
much closer. The probability for chance coincidence between the 
March 5 burst source position and Nl.i9 has recently been estimated®® 
tc be <10"®. 



The synchrotron coolinc niodel removes one of the principal 
difficulties posed by a burst source in the IHC, namely that of the 
very high brightness temperature ("lMeV)of the observed radiation, 

For this model, the bulk of the observed continuum emission, rather 
than being therml radiation of particles with kT i 100 keV, is 
optically thin synchrotron emission of the "MeV electrons and positrons, 
These particles, after losing their kinetic energy by synchrotron 
radiation, annihilate and produce the observed "O.liS MeV linr. The 
minimum magnetic field required for a cooling time shorter than the 
annihilation time is of the order of a few times 10^^ gauss. This 
value turns out, also, to be the minimum B required to confine the 
e'*'-e” plasma that produces the observed annihilation line. This 
confinement is crucial, because otherwise the super-Eddingtonian 
burst luminosity would cause relativistic expansion of the emission 
region in clear conflict with an observed redshifted line. 

An important addition has been recently proposed^** to the synchro- 
tron model for the March ij transient. The MeV e‘*’-e’“ pairs, in addition 
to producing synchrotron radiation, could also Compton scatter their 
own synchrotron photons and thus produce the observed continuum 
above '0.5 MeV. Since both the synchrotron and Compton components 
of the observed spectrum depend linearly on the same relativistic 
pair density, the model. allows a distance determination. The 
distance d depends only on B, the mean relativistic particle energy, 
^mc^ euid the area of the emitting region, A. For B'2xl0^* gauss, 

^mc^ ' 1 !»lev, and A of the order of a neutron star surface area, d 
turns out to be approximately the distance to the LMC. If the source 
were much closer, Compton scattering of the synchr'tron photons could 
not explain the high energy continuum vinless the emitting area were 
very small. For example if d ' 100 pc, A would have to be only a few 
millionths of the neutron star surface. 

In a separate paper in this volume, by Ramaty, McKinley and Jones, 
the possibility of gamma ray amplification through stimulated annihi- 
lation radiation (grasar) has been proposed. This process could 
produce a narrow annihilation line without requiring a low temperature 
annihilation region. Grasar action requires an e'^-e“ pair density 
which substantially exceeds the thermodynamic equilibrium density at 
the temperature of the burst source ( '10^®cra~^at T”3xl0^K). This is 
equivalent to the population inversion of a regular maser. But it is 
not clear yet how such an inversion could be produced in a gamma ray 
burst source. If an inversion can be produced, however, grasar action 
would lead to a narrow emission line at 'O.I 43 MeV without the require- 
ment of a large gravitational redshift. One strong observational argu- 
ment for grasar action would be the detection of narrow annihilation 
lines of width less than abost o.lMeV. 

Nuclear Line Emission 

At least one gamma ray burst, that of November I 9 , 197? also 
shovs^^ evidence of possible nuclear line emission. Its spectrum, 
reviewed by Teagarden in this volume, shows a strong narrow line at 
about 0.7*^ MeV, and other features at higher energies. The continuum 
spectra of gamma ray bursts are suggestive of optically thin 
bremsstrahlung emission of electrons with temperatures ^lO^K. 
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Calculations the expected nuclear deexcitation and 
radiation capture line emissivities from thermonuclear reactions in 
hi^ temperatu:'e plasms have been made for a temperature range of 
10 ® to lo'^K. 



Figure 5. The emissivities of the principal nuclear deexcitation and 
proton radiative capture lines resulting from thermonuclear reactions 
in a plasma of local interstellar composition. ITiese lines result from 
the follovinp reactions: 5.I49I4+ MeV from 2H(p,v)^He; 2.370 MeV 

from y)13n; 1.63U MeV from I4.I438 MeV from 6.129 

MeV from and 0.81*7 MeV from ^^Fe». 
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The calculated emissivities of these lines from a plasma of 
local interstelliu' composition are shown in Figure 5. As can be 
seen, the stron^'est lines for such abundances at tenQieratures of 
10® to ‘3xl0^K are radiative capture lines from ^H(p,v)^He and 
^^C(p,-Y ) ^ ®N. Proton capture on leads to a line at the reaction 
Q' 5.^9** MeV plus the Kinetic energy available in the center of 
mass of the interaction which depends directly on the ion temperature. 
Proton capture on is dominated by a strong, narrow resonance in 
the capture cross section at about 0.li6 MeV/nucleon, leading to a 
line at '2.37 MeV whose energy is relatively independent of the ion 
temperature. 

For -ion temperatures of about -3x10^ to 3x10^ nuclear deexcitation 
lines should be the most intense. The principal lines from a plasma of 
local composition are those at 6.129 MeV from U.U38 MeV from 

‘^C*, 1.63*t MeV from ^®Ne* and 0.81*7 MeV from ^®Fe*. These lines should 
be relatively narrow, thermally broadened to widths FWHM at ~10^®K of 
about lUo, 120, 30 and 10 KeV respectively. 

At ion temperature ^ 3x10^ neutron production exceeds the 
emissivity of the deexcitation lines and neutron radiative capture 
lines could become the most intense. If the capture occurs in 
the hot plasraa(;i 3x10^ ®K) ir which the neutrons are made, however, 
the capture lines would be greatly shifted and broadened by 
the large available energy in the center of mass. In particular 
the r.223 MeV captxire line from Hl(n,y)^H and the 7*632 and 
7.6i*6 MeV lines from ^®Fe(n,y )®"^Fe could be shifted by a few 
MeV or more and broadened to a width of more than an MeV. 

The '0.7^ MeV line and a possible line at '1.1 MeV, reported by 
Teegarden in this volume, in the spectrum of the November 19, 1978 
burst could be nuclear deexcitauion lines at 0.847 MeV and 1.238 
MeV from the first and second levels of ^®Fe, redshifted by z'0.l4. 

Such a redshift is also implied by the broad feature at '0.42 MeV, 
if it is redshifted e'*'-e" annihilation radiation. This redshift is 
within the range of gravitational redshifts expected at the surface of 
neutrons stars . 

But the high relative luminosity of the '0.74 MeV line, " 9 % of 
the total burst luminosity observed above 200 KeV, puts strong 
constraints on a nuclear origin of the line. In particular, if this 
line is redshifted 0.847 MeV emission accompanied by an optically 
thin, thermal bremsstrahlung continuum, then the line-to-eontinuum 
luminosity ratio requires that the emission come from an iron-enriched, 
two- temperature plasma with Fe/H - 10”^, an ion temperature 
> 2x10^ ®K and an electron temperature ^ PxlO^K. Such a plasma of 
vmit optical depth on the surface on a neutron star at a nominal 
distance of 100 pc, emitting '10^® erg see"^ in bremsstrahlung and 
'10^^ erg sec“Mn ^^'Fe deexcitation lines, could produce the observe'd 
burst spectrum. 

\n alternative suggestion"^ that the '0.f4 MeV line is a strongly 
redsh: fted (2'0.3B1 line at 1.022 MeV from single photon annihilation 
of pairs in very strong magnetic fields ('lO^^gauss) has not 

beeri studied in any detail but such a process could not account for a 
possiVle line at "1.1 MeV or higher energy features. 


GAMMA RAY LINE TRANGIENTG 


There is apparently another class of ray transients in 

which essentially all of the observed radiation is in the form of line 
emission. Such a taarama ray line transient was discovered by Jacobson 
et al.^® with a high resolution Ge detector on June 10, 197^ from an 
unknown source. This event, lasting about twenty minutes, was charac- 
terized by strong emission in four relatively narrow energy bands at 

0. 1 4 0-O.I +2 MeV, I.T^^-1.86 MeV, 2.18-2.26 MeV, and 5 . 9 I 4 - 5.96 MeV with 
no detectable continuum. A detailed description of this observation 
is given by Ling et al, elsewhere in this volume. 

As pointed out by Jacobson et al.^® the line identification alone 
appears to require strongly redshifted emission, since the lines at 
'O.Ul MeV and '5.95 MeV cannot be identified with any unshifted lines 
that would not be accompanied by much stronger companion lines which 
were not observed. All of these lines however can be identified with 
lines from the most intense nuclear emission processes expected in 
nature, positron annihilation and neutron capture on hydrogen and 
iron, if we assume that the emission comes from a gravitationally 
redshifted region with a z ranging from '0.2 to 0.205. In particular 
the 5*9^-5.96 MeV line can be identified with the two strongest lines 
from neutron capture on the iron crust of a neutron star having a 
surface gravitational redshift of z - 0 . 285 . Since ®®Fe is expected 
to be the dominant constituent of neutron star crusts this line alone 
strongly suggests a neutron star source for the transient emission on 
June 10 , 19Tii . 

The O.UO-0.1^2 MeV and the l.T^-1.86 MeV lines can similarly be 
identified as the positron annihilation and neutron capture in a 
hydrogen atmosphere of density nj^ > 10 ^® Hera”® extending from the sur- 
face z of 0.285 up to an altitude corresponding to a z of about 0 . 2 . 
Such an atmosphere is presumably a temporary feature perhaps formed 
by infall of gas from an accretion disk. This minimum hydrogen density 

1. s required in order that a significant fraction of the neutrons are 
captured before they decay. The width of the 'O.^i MeV line also places 
a limit on the average temperature of this atmosphere of <2 x 10 ® K. 

The unshifted (z 1 Oj^.O?) line at 2.22 MeV requires a comparable high 
density region well above the neutron star surface and can be under- 
stood as neutron capture in the atmosphere of a binary companion. 

A surface redshift of '0.205 implies a neutron star mass of about 
I.I 4 to 1.8 M,and a radius of 10 to 13 km, depending on the assumed 
equation of state. The observed 25keVfull width of the '5*95 MeV line 
reduced by the instrumental broadening of 5 keV and the redshifted ®®Fe 
doublet separation of 10 keV gives em emitted width of lOkeV. If the 
line were produced unformly over the surface of the neutron star, this 
width would give a lower limit for the rotation period 

T, _ Fy i4TTr sine sin6 

^ ’ Je J ' T" 

where 0 is the angle between the rotation axis and the direction of 
c bservat ion . 
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If the lines do result from neutron capture near the surface of a 
neutron star, however, the observed line fluences would then require 
either rapid mixing or rapid breakup of the resultant and ^'Fe formed 
by the capture since these nuclei would quickly accumulate to many 
optical depths for any reasonable source distance. The column density 
of and ^^Fe, produced at the surface of a neutron star during the 
June 10, 197*4 transient, would amount to ‘’3xlO^°D^OO 
'lO^^D^OO ^^Fe cm*"^ where I^ioo assiuned distance in units of 

100 pc . 

But even assuming such an identification for the lines, the origin 
of the neutrons and positrons required to produce them is still 
problematical. Although one might imagine that it should not be hard 
to find ways of generating neutrons directly on the surface of a 
neutron star, the observation of the unshifted neutron capture line 
on hydrogen would require that the average energy of the neutrons 
produced at the neutron star surface must be '’150 MeV, equal to the 
escape energy from the star, so that a significant fraction could 
escape to a lower z region to be captured. There is no obvious process 
for generating such energetic neutrons. It also seems unlikely that the 
neutrons and positrons could have been produced in nuclear reactions 
in a flare on the binary companion because of energetic consideratior.i . 
One possibility, however, would appear^® to be that the neutrons and 
positrons which were captured and annihilated in the vicinity of the 
neutron star and its binary companion were all produced by nuclear 
reactions occurring when gas, flowing episodically from the companion, 
was gravitationally accelerated to n energy of > 1 MeV and collided 
with gas in the outer part of the accretion disk around the neutron 
star. In addition to the observed lines we would also expect other 
lines and continuum “rrission, but the expected intensities are not 
necessarily inconsistent with observations. 

Lastly, since the sum of the neutron-capture gamma-ray line 
intensities for the transient event of June 10 , 197 ** was ~ 3x10” 
erg cnr^ sec”*, the neutron capture line luminosity of the source 
i sec~^ . Assuming a maximum neutron yield of * 

< 2x10^ neutrons erg”'^ , see ref. (36), and an average neutron capture 
photon energy of ‘ 3 MeV, the total source luminosity L must have 
been ^ lO^^D^j^p erg sec”*. 
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A REVIEW OF THE 1979 MARCH 5 TRANSIENT 


T. L. Cline 
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NASA/Goddard Space Flioht Center, Greenbelt, MD 20771 


ABSTRACT 


The understanding of the 1979 March 5 event remains a problem 
of central Importance for researchers In gamma ray transient 
astronony. Efforts In the study of the observational results and In 
the development of Interpretations regarding the nature of this 
event and Its relationship to the variety of other gamma ray burst 
phenomena have continued through the past 2 years. A consensus of 
opinion has not yet been reached regarding Its possible origin In 
N49 at 55 kpc distance, versus In an Invisible source 3 or 4 orders 
of magnitude closer, although Interpretations favoring N49 appear to 
be presently gaining momentum. This presentation outlines the 
existing data In a review of what remains the most singular high- 
energy astrophysical phenomenon of the spuce age. 

INTRODUCTION 

1 assume that most of the persons Interested In this subject 
are familiar with the March 5, 1979 gamma ray transient; since the 
bulk of the observational data has been published for some time w;v;» 
review seems at first unnecessary. However, progress has been and 
continues to be made In both the experimental and the Interpretative 
reanalyses of the measurements. Thus, although the event remains an 
historic fact. Its study Is very much a current activity; In fact, 
the ‘correct* Interpretation Is undoubtedly yet ahead of us. I 
personally have found this event to be fascinating, both because of 
Its apparent observational rarity and because of Its role as a 
generator of theoretical Ideas and related calculations. Will we 
see a gamma ray transient like this again In our research 
lifetimes? Probably not, at least, probably not when we shall be as 
prepared with third-generation Instrumentation as we happened to be 
with second generation Instrumentation on March 5th of 1979. 

However, If detector sensitivity and resolution can be sufficiently 
Improved, March 5-1 Ike events of much weaker visual magnitude may be 
observable and capable of detailed study, v/hether arriving from more 
distant stellar or from extragalactic sources. Will an entirely new 
modelling (such as the grasar, to be unveiled at this conference^) 
be found to be applicable? If so, some new exper1nH>nt£l approach 
Involving the measurement of, for example, polarization, or 
submillisecond oscillations, or gairnia ray coherence, may provide 
another observational break though. In addition to that which gamma 
ray line spectroscopy now provides. 

We can first review the question of whether the 1979 March 5 
event Is so different from other gamma ray bursts as to be termea 
'unique*, which I have described It to be In an earlier review.^ 

Invited Review Paper presented at the La Jolla Institute Workshop 
on Gamma Ray Transients, 1981 August 5-8. 
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The variety of gamnna rty burst data we have been shown at this 
conference. In fact, leads one to believe that every gamma rty burst 
Is unique! However, the first basic difference chat appears to 
confront us Is that while the typical or so-called 'classical' gamma 
ray bursts seem generally to have time histories of random 
structure, with randomly evolving spectral features, the March 5 
event appears to be especially ordered, with a singularly fast rise 
time, a single very Intense peak, and a subsequent, well-defined 



Fig. 1. Time 
history of the 
1979 March 5 
transient as 
observed with 
the I SEE-3 
space probe^. 


oscillation. The other basic difference Is the fact that the March 
5 event Is the only one with a candidate source object. Of course. 
Its features may only seem to be atypical due to their greater 
resolution afforded by the great apparent brightness of the event; 
also. Its source Identification could be a cosmic accident. To 
'write off the event In this manner would be to completely Ignore a 
possibility of great value to astrophysics. In ity opinion. 

The March 5 event possesses a 420-keV spectral feature^ In 
common with some other classical bursts^*”; this, together with the 
regular oscillation (which Is not clearly exhibited by classical 
events) has prompted the general belief that both originate In 
netitron star phenomena. However, they may of course originate by 
means of two or more entirely differing mechanisms. In fact, the 
420-keV feature also has a novel and entirely differing Implication: 
as we shall hear later In this conference^, the neutron star 
redshift explanation Is no lon^r the only possible Interpretation 
of the 420 keV feature. Thus either the March 5 or the classical 
event connections to neutron star redshifts or both may yet be found 
to be Incorrect or Imprecise. 

It Is a fair statement, unfortunate though It may he, that no 
mathematical treatment has yet been evolved to < ystematically 
classify or categorize gamma ray bursts; thus, the probability of 
uniqueness of the March 5 event has not been quantitatively 
estimated. To pursue this, we can look at some of Its features In 
greater detail. 


RATE (COUNTS /SEC) 
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OSCILLATIONS 


Of the several unusual features of the March 5 event* Including 
the rise* the Intensity* the oscillations* the spectrum and the 
direction* the oscillations were Invoked first as providing 
basically new Information/ Except for a weak but suggestive 'o 




No. i 

No. 2 

No. 3-4 
No. 5-7 
No. 8-11 
No. 12-16 

No. 17-22 
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Fig, 2. (Above) Venera-11 observations 
of the March 5 event (horizontal scale 
5 seconds/mark), Indicating the first 
clear evidence for periodicity In a gamma 
ray transient (from Mazets et al.^). 


Fig. 3. (Left) ISEE-3 observations of 
the first 22 cycles of the March 5 event, 
plotted on an 8-second per period basis 
folded with an Increasing number of cycles 
per plot, to compensate for the decreasing 
statistical validity per cycle. The large 
pulse and the secondat 7 features remain 
constant In time, yielding an average 
period of 8.00 t .05-second*^ 


4-second, - 5 cycle long, periodic effect In the time history of the 
29 Octt^r 1977 gamma ray burst event,” no evidence Is known for apy 
cyclic features in classical gamma ray bursts. The clear 8-second 
periodicity In the March 5 event provided the first evidence of that 
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kind for the neutron star origin model of a gamma r<^ transient. 
(Spectral evidence existed before that time for one very si ok gamma 
ray transient^ although It was not then and still Is not clear 
what relationship that Isolated phenomenon has to classical gammr 
ray bursts. Of course. It also Is not clear what relationship the 
March 5 event has to classical events.) Figure 2 shows the results 
from the spaceprobe Venera-11 as published by Mazets and his 
colleagues of the Leningrad group*^*' Both the general 8-second 
repetition and the ^ 4-second modulation are clear. Similar time 
histories were exhibited by a companion sensor on Yenera-12. The 
oscillations were tracked for about 3 minutes with the Goddard 
experiment on ISEE-3, as seen In the time history In Figure 1. but 
of course suffered from decreasing statistical significance as a 


Fig. 4. The time history 
as observed with the 
Franco-Soviet Instrument on 
Venera-12. Indicating 
evidence for 1-second 
features. *2 


i m U 

function of time. This effect Is compensated for in the presen- 
tation of Figure 3. In which an Increasing number of cycles are 
Included per plot. Her« the constant phase of the comoound 
structure Is evident, providing a period of 8.00 ± 0.05 second.^ 
This periodicity Is clearly not a fluctuation but very definite 
evidence for a cyclic phenomenon; whether the 8-second feature Is 
rotational, processional or radial Is another question. The fact 
that the main and secondary pulses decay differently has been 
Invoked as evidence for a directional emission from a rotating 
neutron star.“ 



Weak evidence also exists for features of time structure finer 
than the 4-second Interpulse. Figure 4 shows the Venera-12 data as 
published by Vedrenne and his colleagues of Toulouse and Estulln and 
his colleagues of Moscow. Here, as In earlier figures, 
suggestive fluctuations are evident on a shorter time scale. The 
spectral power as a function of period derived from these results. 
shown In Figure 5. and the power spectrum as a function of frequency 
of 3 minutes of Ploneer-Venus Orbiter data^^ of Evans and his 
colleagues at Los Alamos, shown In Figure 6. both Indicate second- 
order structure In the - 0.7 to 1.1 -second region. (The I SEE-3 
results do show some detail but Its time history may be somewhat 
spin-modulated, although not severely, since the source direction Is 
only 3.5 degrees from the ISEE-3 spin axis.) Finally, one can 
suspect that even finer time variations could exist within the first 
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oscillation. In Figure 7 the Leningrad Veneri*ll data hint at a - 
0.2-second fluctuation. 




Fig. 5. Spectral power 
evidenced In the Franco- 
Soviet Venera-12 results. “ 


Fig. 6. Power spectrum from 
the Ploneer-Veous-Orblter 
measurements.*^ 





Fig. 7. Detailed Venera-11 
time history of the early 
part of the event, 

Including the first Intense 
peak, Its decay and the 
buildup of tbe first 
oscillation.^*’ 



6 


BURST TIME HISTORY 

The Initial hi gh> Intensity burst of radiation In Itself may be 
even more anomalous than the oscillating portion of the March 5 
event. Its features, not shared by the general population of 
classical gaMM ray bursts, are a fast rise, a single peak, a smooth 
decay and an Instantaneous Intensity near detector-saturation for 
most of the Instruments (of the order of one count per 10 
microseconds). In fact, the unknown spectrum of hard X-rays below 
detector threshold for these gaimna ray detectors presents some 
uncertalnity of possible measurement distortion; It Is hoped that 
the effect Is not major, a belief supported by the fac . that the 
observations generally mutually agree. For exa^le. the onset rise 
Is faster than can be measured with every Instrument, yet the peak 
counting rate Is seen to occur usually about 20 milliseconds Into 
the event. (Unfortunately the ISEE-3 y-ray spectrometer falleo 2 
months earlier, depriving us of a comparison measurement made with 
germanliMi rather than with scintillators.) This onset consists of a 
sudden Increase of several orders of magnitude within the resolution 
time of. for example. I millisecond In the ISEE-3 Goddard/1^1 
scintillator. 3 as shown In Fl^re 8. The exponential time constant 
of Intensity increase Is therefore less than 200 microseconds. 
Implying a light travel distance of less than 60 km. 

Fig. 8. (a) Onset of the high- Intensity 

portion of the 1979 March 5 transient.^ A time 
constant of less than 0.2 ms Is Inferred from 
the Increase of two orders of magnitude from 
near background to essentially full Intensity 
within a resolution time of 1 ms. (In this 
Instrument the time to accumulate 64 photons Is 
recorded to I m accuracy; the first several 
readings are In fact 1 ms and 2 ms accumu- 
lations.) This < 1 ms full rise In the onset 
shape Is^alsg seen with Ploneer-Yenus- 
0rb1ter.3. » 

No Other gamma rav transient exhibits as brief a rise time, within a 
factor of 10 to 30 depending on the various detector limitations 
consldereo. (For example, the Goddard Helios-2 detector Is 
Instrumented with slower. 4 millisecond temporal resolution* the 
entirety of the fast spike as observed with that Instrument^ Is 
shown In Figure 9. Four years of gamma ray bursts were monitored 
before the Hel1os-2 spacecraft, was turned off; only the March 5 
event produced a tr1gg<er response In the 4-ms circuits; all other 
events triggered In Its 32 or 2SO-ms modes.) However, other events 
have been ^served that are briefer than the 150-m1111second total 
duration observed here. Two events of moderate Intensities were 
detected In 1979 with Venera and PVO. one of which was also observed 
with Venera-12 giving a directional determination at high galactic 
and celestial latitude (k. Hurley. R. Klebesadel. private 
communication). The question of whether such fast, single-spike 
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events occasionally detected* Including perhaps 3 more In the 
archived Vela data (R. Klebesadel, private communication)* are 
merely the tips of more complex classical-event Icebergs or are 
March 5-1 Ue events* cannot yet be answered. 



Fig. 9. The 150-m1111second wide 
Initial burst. A transition from 
a slow decay of ~ 150 ms time 
constant to a steeper decay of 
35 ms time constant Is seen about 
100 ms Into the event. These data 
are plotted on a ~ 4 ms per 
accumulation time basis. ^ 



Fig. 10. The data* on a ~ 2 ms 
per accumulation basis* from the 
Franco-Soviet experiment on 
Venera-11* collected during the 
Initial one-quarter second of the 
March 5 event. 


The question of very fine time variations In the March 5 event 
time history Is Interesting for several reasons. (These range from 
the existence of 0.1 to 100 ms oscillations, due either to an 
Inherently fast neutron star spin period or to model -dependent 
transition damping vibrations,*® to the possible existence of a very 
short duration nonrandorolcity In the counting rate that would relate 
to a possible coherence, as In certain other models. M Detailed 
studies of the first 125 milliseconds of the event, examples of 
which are shown in Figures 9 and 10* have been carried out using 
data from all available high time resolution sources* Including the 
Einstein Observatory X-ray monitor. The Venera time histories hint 
at a 25 msec effect with marginal statistical accuracy (K. Hurley* 
private communication) and the Einstein time history* of only 55 
msec duration* can be Interpreted to be consistent with two 27 msec 
features (M. Welsskopf* private communication); these departures 
from randomicity appear to be tantalizing but Inconclusive. 



8 


INTENSITY AND SPECTRUM 

The total Intensity of this burst Is that of a typical, strong 
gamma ray burst, but the Instantaneous inU'nsity of the peak is by 
far the greatest yet observed. Figure 11 shows a scatter plot of 
event total intensity versus maximum instantaneous intensity> using 
the sample of Helios-2 events detected from January 1976 to the 

Fig. 11. Diagram of total 
versus maximum intensity for 
the gamma- ray bursts observed 
with Helios-2 from 1976 up to 
the 1979 March 5 event. The 
more extensive burst observa- 
tions of the Vela system 
cannot be used for this 
purpose since they employ a 
geometrically lengthening 
time base whereas bursts do 
not always have maximum 
intensity near onset. If a 
search is made for Vela 
transients of March 5-like 
character, however, a 
comparison of total versus 
onset intensity gives a 
similar result. 

March 5 event. It is also the case that this is the most intense 
event monitored with the Vela satellites for the decade beginning in 
1969 (R. Klebesadel, private communication). The observed intensity 
is at least several times 10"^ erg cm*‘ sec'*, a lower limit due to 
the unknown fluxes below the 30 to 50-keV instrument thresholds and 
to the unknown effects of pulse pile-up at these energies {consi- 
dering that most instruments were running near saturation). The 
intensity of the oscillating portion is - 10“^ that of the peak; 
since this event is so intense relative to most gamma ray bursts, it 
is possible that some of the classical events of average or weak 
intensity could possess an analogous but unobservable periodic decay 
con^onent. 

The spectrum of the March 5 event is another of its 
anomalies. Although both this event and some classical events have 
the famous 420-keV featurei^.5.6 this event has a considerably 
softer spectrum below ~ 100 keV than do typical classical bursts. 
Figure 12 illustrates the spectra of the peak of the March 5 event 
and of the oscillating portion,^ compared with that of a typical 
classical event. The 420-keV feature is clearly evident in the 
intensity peak (the loss of the ISEE-3 germanium spectrometer, 2 
months earlier, deprived us of a knowledge of the inherent width of 
this peak). Besides the excess in the 30 keV region and the 420-keV 
feature, the other well-known aspect of this spectrum is the lack of 
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Fig. 13. Ploneer-Venus-Orbiter measurements of the spectral hard- 
ness of the March 5 event (1n the above/below 100-keV region), 
Illustrating variable behavior at the onset and shoulder of the 
Intensity spike, and ■ in phase with the first two oscillations. 


higher energy photons, particularly above 1 MeV. This deficiency is 
probably due to the fact that photons above the pair production 
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threshold energy would be, of necessity, removed by that process due 
to the extreme density at the source, particularly If the source Is 
very distant. 1^ The spectrum of the oscillating portion of the 
event Is even softer than of the onset, with no detectable counts as 
high In energy as 400 keV. The only other published spectral 
observations of this event are shown In Figure 13, Illustrating a 
spectral hardness factor (counts > 100 keV/50-100 keV) as a function 
of time. 11 The hardness Is maximum at the onset and at the time of 
the change of rate of Intensity decay, or shoulder, at - 100 msec 
into the event. As shown In the second Illustration, the spectral 
hardness In the periodic portion Is less than In the onset. In 
agreement with the Leningrad results, but Is shown to vary In phase 
with the 8-second oscillation, with higher enerw photons in greater 
quantity during the times of greater Intensity*^ 

RECURRENT EVENTS 


The March 5 event was the first gamma ray transient to exhibit 
recurrent events, I.e., to be associated with repeated bursts having 
a common source direction. Three events of very weak Intensity, 
detected only with the Leningrad experiments on Venera-11 and -12, 
were observed to follow the March 5 event by delays of 0.60, 29 
and 50 days with Intensities - 3, 1 and 0.5 percent that of the 
March 5 Intensity, respectively. 1^ Their time histories, shown in 
Figure 14, are not the same (on a reduced scale) as the March 5 
profile but generally slower, although their spectra, as seen in 
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Fig. 14. Counting rate histories of 
the delayed March 6, April 4 and 
April 24 events, found to have source 
directions consistent with that of 
the March 5 event. 15 


Fig. 15. Spectra of the 
events of the previous 
figure, seen to be not 
unlike that of the decay 
portion of the March 5 
event. 15 
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Figure 15, are similar to that of the softer, oscillating portion 
although statistically limited due to their very low Intensity. 

Their olrectlons agree with that of the March 5 event (the Venera-11 
to Venera-1? time delays provide an accurate one-dimensional source 
measurement In conjunction with a rough source direction estimated 
from the directional count rate characteristics). Given their 
sequential connection to each other and to the March 5 event. It Is 
entirely reasonable to assume this to be the discovery of a common 
burst source emitter. Only one other series having a common source 
direction has been detected;^' however, that 3-event series does not 
follow any intense primary or 'parent' gamma ray transient. Its 
source direction is In the galactic plane at - 45 degrees galactic 
longitude, consistent with either a nearby or a very distant 
galactic source object. 


DIRECTION 

One remarkably fortuitous historic accident was that of the 
occurrence of the March 5 event during the October 1978 to December 
1979 time frame. This is the interval when the interplanetary gamma 
ray burst network possessed its full instrument complement and was 
positiorod in trajectories far from the Earth for their maximum 
Interspacecraft event timing resolution. The identification of the 
supernova remnant N49 in the neighboring galaxy of the Large 
Magellanic Cloud as a possible source was made as one of the 

first accomplishments of that network. This measurement made 
possible the enlarged scientific controversy over the nature of the 
physical process involved by immediately creating the distinct 
possibility that the source object was at a 55-kpc distance, two to 
three orders of magnitude farther than would have been assumed had 
the directional information consisted instead of the formerly 
available resolution. That would have been several square degrees, 
providing only rough overlap of the LMC and various stellar regions 
In the constellation of Dorado at 30-odd degrees negative galactic 
latitude. Figure 16 illustrates the source location, and shows the 
extent of the densest optical part of N49. Another aspect of this 
Identification is that, by contrast, high-precision directional 
studies of classical gan»na ray bursts made with this network have 
provided source f1g]ds that are optically empty down to at least the 
18th magnitude. An additional factor is the fact that N49 is 
an X-ray emitter, strengthening the case for its Identification, 
since one might assume there would be some X-ray /gamma- ray corre- 
lation. Finally, it is also the case that not only other precise 
burst source locations but previous larger source field determi- 
nations also found no correlation with either transient or steady 
X-ray sources or other objects that could be obvious source 
candidates'^, 23, Figure 17 shows the source position plotted on the 
X-ray contours of N49 measured (quite accidentallv) within a few 
days of the event with the Einstein Observatory.^^ A remarkable 
aspect of this X-ray survey, however, is the fact that no point 
source was found: this does not necessary provide a weakening of 

the association of the event with N49 since It yields an upper limit 
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(at 10-9) to the ratio of point source X-ray strength to byrst 
strength independent of source distance. (The X-ray measurements 
also Included a comparison, made by chance before and after the 



Fig. 16. Source location of the 
March 5 event shown on a reproduction 
of optical field, In which the 
densest portion of the N49 snr 
nebula is visible. 


Fig. 17. Source location of the 
eventl9 plotted on the X-ray 
surface brightness contour map of the 
N49 and (N49) region, as observed 
with the Einstein Observatory high- 
resolution imager. 24 Here the 
nebula is seen as 2 arc minutes wide, 
and the nearest neighboring nebula, 
(N49), is resolved at the weakest 
intensity level. No point X-ray 
source was resolved. No chanae in 
X-ray intensity above - 2x10"^^ 
erg cm-2 s-1 was observed from 
shortly before to several days after 
the event. A comparison X-ray study 
with Einstein has been recently made 
to search for a slow, delayed effect, 
but the results are not available at 
this time. 


burst, showing again an upper limit of - 10-9 to the intensity 
change ratio.) Thus, the fact that N49 is an extended object, 
suggesting a greater likelihood of chance association, is contrasted 
by the facts, supporting its identification, that it is a supernova 
remnant - an object associated with neutron star phenomena, and that 
it is an X-ray emitter. Further, the visible stars in the source 
field have been determined to be also at the same distance as the 
LMC.25 thereby removing their distraction from this association. Of 
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course, the argument for chance association Is only a factor of two 
larger than if the candidate object were a point object, since the 
error box and N49 consist of about the same solid angle. That 
probability, considering the total solid angle of extended X-ray 
emitters and the total number of point X-ray sources (times the 
source box size) Is between lO"^ and 10’®, too small to be a strong 
argument against the Identification. The basic aroumen'. against N49 
Is of course Its distance and the 10^5 sec’^ emission iruplled 
If the radiation In this event Is Isotropic. 

DISCUSSION 

The observations of the March 5 event, although Inconclusive 
regarding the question as to whether this event Is so different 
(that It must have been produced by an entirely different 
mechanism), clearly provides support for the hypothesis that Its 
nature and possible origin at N49 are phenomena so unusual that they 
deserve critical study per se (rendering the question as to whether 
Its understanding will enl ighten us as to the nature of ‘classical* 
gamma ray bursts, or vice versa, as secondary), A variety of early 
papers concerning this event tended to dismiss the N49 association 
as chance, avoiding the very real problems In the theoretical 
contortions necessary to get around the photon self-absorption 
Inherent In fitting a physically possible mechanism to 
sec*l luminosity required by an Isotropic N49 emission. 

It had been known for some time that gamma ray burst Intensities 
Implied a several hundred pc source distance If treated at face 
value, given that the spectra entended beyond the pair production 
threshold;16.29 recent updates show that, given nearby classical 
burst source d1 stances, 3^ the March 5 event source distance would be 
accordingly only -2 pc, keeping the optically thin Bremsstrahlung 
treatment.31 Creative treatments of the March 5 Issue, taking the 
N49 association as a serious and possibly fruitful possibility, have 
only recently started to gain momentum. 

The first detailed calculations of the March 5 event spectrum, 
treating the process of pair production, annihilation and scattering 
with an optically thin synchrotron mechanism (In the Intense 
magnetic field of a neutron star requiring a - 20 percent gravita- 
tional redshift), fit the observations surprisingly well. Figure 18 
shows the calculations of Ramaty and his coworkers, as outlined In a 
November 1979 gamma ray burst conference. 32 jh^ low energy excess 
emission, the position and width of the line anu the general shape 
of the continuum are all seen to be duplicated. In essence. The 
more worrisome problem of the source mechanism was first addressed 
In a sequel to this work. In which the same group showed that an 
earlier model of the storage of energy In neutron star vibrations. 
Involving the release of (undetectable) gravitational radiation, 
provided a good fit of the vibrational dandling time to the 150 msec 
width of the Initial burst (see Fl^re 19). This model provided the 
necessary storage mechanism that made the 55-kpc distance physically 
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Fig. 18. The calculated spectrum 
of the intense portion of the 
March 5 event, based on an 
annihilation-creation model 
Incorporating synchrotron losses 
and a redshift effect In a 
neutron star environment, 32 
compared with the observations. 4 
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Fig. 19. The quadrupole 
gravitational radiation 
damping time versus gravi- 
tational redshift for neutron 
stars, with the dashed line 
connecting calculated values 
having the same equation of 
state, compared with the 
March 5 event. 


possible, and linked the Internal, gravitational process to the 
external, magnetospheric, observable phenomenon. More recently, the 
Rama^ model has been extended by Liang to Include the higher energy 
effects of the Inverse Compton process. 33 These calculations not 
only fit the observed spectrum a little better, but also provide a 
derivation, from first principles, of the luminosity at ~ 1044 erg 
sec'l, giving an order-of-magnitude fit to the N49 source. 1 do not 
Illustrate It here since It is presented in this Conference. 33 
Given that these calculatlonal exercises provide consistency, or 
better, with the N49 source association, we are freed from the 
necessity of assuming that It must be a chance association. Final 
proof of which view may be correct Is quite possibly a long way off, 
however. 


An ongoing reanalysis of the March 5 event directional data has 
provided a more precise source location; 34 this new error box Is 
Inside M49 although not at Its center (see Figure 20). The 0.1 arc 
m1n2 size of this source field Is only about § percent that of the 
original conservative treatment of essentially the same data, making 
searches for point objects possibly a lot easier, but reducl.ig the 
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Fig. 20. The revised 1979 
March 5 event source location, 
entirely consistent with N49 
but In an eccentric location, 
15 to 55 arc seconds from Its 
center. 3^ Given an unknown 
radial component, a possible 
location of the source may 
even be at the edge of the 
shell. This Is the most 
precisely determined gamma 
ray source location In 
existence. 


probability of chance identification with N49 by only a factor of 
two (since the N49 solid angle remains 2 arc min2). This result 
would imply a motion of the neutron star from the center of the 
remnant, assuming N49 to be source, of - 900 ± 400 km sec"% if the 
age of the remnant is 10,000 years. This is probably not an 
entirely unreasonable value. Of course the neutron star involved in 
the March 5 emission may not be the parent neutron star of the N49 
snr itself, given the snr and neutron star density in the LMC. 

What then is the relationship of this event to the bulk of 
gamma ray bursts - why Is there a gap in the size spectrum? No one 
would extend the N49 association to the point of attributing the 
March 5 event, as the brightest of all bursts, to an LMC origin 
because that galaxy is the nearest of the external galaxies that, in 
turn, produce the weaker. Isotropically distributed, classical 
bursts. This size spectrum difficulty is avoided in one way by 
assuming that March 5-like events are rare, either intrinsically 
rare, as are supernovae themselves (perhaps for reasons that are 
related), or apparently rare, because of a beamed emission 
anisotropy (or a mixture of both); in this case the detection of one 
event in a decade from the neighboring galajy LMC is not really less 
likely than from our own galaxy. 35 (The nearby origin model has no 
corresponding way out of the size spectrum difficulty, other than to 
invoke a simple accident of proximity, like the accident of 
coincidence with the direction of N49 itself.) Continuing this 
speculative notion, a possible link between the March 5 event and 
the other bursts is the series of three weak, delayed events, 
compared with the existence of another series of three small common- 
origin events that does not itself follow any identifiable large 
event. The source direction of that series, as noted above, is in 
the galactic disk at - 45 degrees galactic longitude. 17 That 
direction is consistent with a galactic source distance of up to 
perhaps 20 kpc, half the distance to the LMC, putting those events 
into a model of consistency with an N49 origin for the delayed March 
6--Apr11 series, except that any primary or ‘parent' event In this 
case was Invisible, due perhaps to an anisotropy of emission. Are, 
in fact, the larger ‘classical' gamma r^y bursts also 'secondary. 
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delayed' emissions, but with repetition times outside the presently 
measurable several -year extents, perhaps linked to their greater 
Intensity* relative to the March 6— April series? This would Imply 
the existence of many apparently Invisible March 5-1 Ike events, 

I.e., of the requirement of a narrow solid angle of emission, 
perhaps several hundred square degrees. Yet some link surely 
exists, defining the occurrence rate of the various kinds of 
transient events from each neutron star source. Finally, one can 
question the basic assumption of the neutron star gamma ray 
transient source concept, accounting for the 420 keV lines as 
evidencing a - 20 percent redshift, gravitationally required for a 
1.4 Mq neutron star, of the 511-keV line. The stimulated emission 
concept Is shown In this conference to produce a 420-keV feature 
without the presence of a gravitational redshift. ^ Clearly a great 
deal of research Is needed to Investigate the emission charac- 
teristics Implied, such as the amount of coherent beaming, that may 
clarify, eliminate or unify these Ideas In terms of the relationship 
to gamma ray phenomenology of the 1979 March 5 event. 


I wish to point out that I have not treated some very recently 
published theoretical studles.pf thegMarcb^B event, but I do wish to 
refer them to your attentlon.'^®*^'*^®*^’*^^ 
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ABSTRACT 


A review of recent results in gaunma-ray burst spectroscopy is 
given. Particular attention is paid to the recent discovery of 
emission and absorption features in the burst spectra^^' ^ ^ 

These lines represent the strongest evidence to date that gamma-ray 
bursts originate on or near neutron stars. Line parameters give 
information on the temperature, magnetic field and possibly the 
gravitational potential of the neutron star. The behavior of the 
continuum spectrum is also discussed. A remarkably good fit to 
nearly all bursts is obtained with a thermal-bremsstrahlung-like 
continuum. Significant evolution is observed of both the continuum 
and line features within most events. 

INTRODUCTION 

This paper is intended as a review of the measurements to date 
of the spectral characteristics of gi»jnma-ray bursts (GRH). It will 
also cover to a certain extent the theoretical interpretation of 
these results where it relates in a direct way to the understanding 
of the observations. Within recent years GRS spectral measurements 
have begun to reveal a very rich phenomenology and to make a major 
contribution to the understanding of the origin and physics of 
GRE's. 

Early measurements^"® were made with instruments designed for 
some other purpose than the detection of gamma-ray bursts, with the 
launching of instruments specifically designed to study the spectral 
characteristics of GRB's^"^^ a great variety of new observations 
have significantly extended our knowledge of the spectral behavior 
of GRB's. In particular the work of Mazets and his colleagues at 
the A F. Ioffe Physical-Technical Institute^" ^ ^ , Leningrad, 
represents the most comprehensive study of GRB spectra presently 
available. Tneir data have among other things revealed the repeated 
- esence of absorption and emission lines in GRB spectra. "Cyclotron" 
sbsorption features are very probably present in a major fraction 3f 
the events measured^*^. In addition, lines at 400-450 keV are seen 
in a smaller number of events that are probably from redshifted 
annihilation radiation''®. Taken together these results are strong, 
if not compelling, evidence for the neutron star oriqin of GRB's. 

In this paper the ea’'ly spectral results are first briefly 
reviewed. Tnen the behavior of the GRB continuum spectrum, its 
form, variability, and evolution, is treated. Finally the absorption 



and emission features, their relationship to other spectral 
parameters and physical interpretation are discussed. 

HISTORICAL OVERVIP-? 

Early results on GRB spectra ^ ^ generally suffered from a 
number of problems; e.g., poor statistics due to small detector 
area, inadequate time resolution and variable intensity modulation 
due to spacecraft rotation. The first published GRB spectra^ were 
derived from a small instrument on the Imp-6 spacecraft intended to 
measure solar flare X-rays and interplanetary positrons. Spectra 
from six events were reported covering the energy range from 0.1-1. 2 
MeV. Tne spectral form for all events was consistent with a simple 
exponential form exp{-E/Eo', with E© " keV. In a subsequent 
paper^ using data from the Imp-7 spacecraft they presented furthei 
evidence that event-integrated GRB spectra could be represented by a 
single functional form. This time the ori'.ial 150 keV exponential 
was modified by the addition of a power-law high energy tail 
(E"2.5), Imp-7 spectra of Cline and Desai^ are reproduced in 

Figure 1. Othei contemporaneous results^ ' ^ ^ ^ ^ were consistent 
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Figure 1 . GK3 spectra from the Imp-7 experiment of Cline and co- 
workers'-. A single functional form, consisting of a ISO keV 
exponential wit>i a power-law tail apnears to fit all the event.s. 
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with this picture. We shall see later, however, that this simple 
picture was, in fact, not correct and that spectral variability from 
event-to-event and evolution within individual events are, in fact, 
typical characteristics of GRB's. 


Of the early GRB spectral results, the most precisely 
determined spectrum was that of the event of 1972 April 27 as 
measured by the y- and X-ray spectrometers on Apollo IG^* This 
spectrum is shown in Figure 2. A combination of two detectors, (1) 
an X-ray proportional counter i.nd (2) an Nal y-ray spectrometer, 
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Figure 2 . Spectrum 
of the Apollo 16 
GRB^^ averaged over 
.he entire event. 


gave the wides*- "nergy coverage (2 keV’ to 5 MeV) of any GRB t>'.us far 
recorded. 'rC'. unfortunate gap between the coveraae of the two 
instruments in the 10-60 keV’ range may have precluded them from 
mahinc the first discovery of cyclotron liner in GRP ' s . In a 
subsequent paper Gilman et al .'*^ showed that the Apollo 1R data 
could be fit remarkably well by a thermal-bremsstrahlunq spectru^i 
with a temperature of 500 keV. It should be noted that this fit to 
the data was quite good over a range of thiee decades in energy. 

The auth.ors naturally assumed that the burst source was, in fact, a 
hot optically thin plasma. TP maintain the observed spectral shape, 
the source must be optically thin to Compton scatterina. In 
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addition, it is possible to place a limit on the size of the 
emitting region the Apollo 16 event based on the most rapid time 
variations observed during the event. Variations as fast as 110 
msec werr “"een, which imply a source size ~ 3000 km. These 
assumptions taken together allow one to place a limit on the 
distance to th*/ source of ~ 50 pc. This early result did, in fact, 
establish a corr^^ct spectral form for GRB’s (i.e., a thermal- 
‘bremsstrahlung-like spectram) . The interpretation, however, in the 
light of more recent data^®'^^ is almost certainly not correct. 

Later results will show that the emitting region is probably not 
optically thin and '’ominrted by free-free emission. 

CONTimiUM SPECTRU M 

The KONUS experiment of E. P. Mazets and coworkers at the A. F. 
Ioffe Institute, Leningrad has yielded by far the most comprehensive 
set of data on the spectra of gamma-ray bursts^' ^ ® ^ ^ • Their 
relatively simple instrumentation consists of six Nal scintillators 
each of ~ 50 cm^ area arranged to provide complete sky coverage. 

The spectrum is measured over the 30 keV-1 MeV interval in 16 
channelr. Identical instruments were carried on board the Venera-11 
and -12 spacecraft. The fortunate combination of all-sky coverage, 
jsence of earth occultation and magnetospheric background effects 
and a 1 1/2-year active lifetime has produced a data sample of ~ 150 
GRB's. The spectrum of the 1970 April 19 GRB as measured by the 
KONUS experiment^® is reproduced in Figure 3. As was true with the 


Figure 3 . Spectrum 
of the 1979 April 
19 GRB from the 
KONUS experiment at 
the maximum of the 
event . 
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Apollo 16 event, this GRB is well fit by a thermal bremsstrahlunq- 
like spectrum. In fact, the radiation temperature of 470 keV is 
close to that of Gilman et al . However, it should be pointed out 
that Gilman et al . incorporated a Gaunt factor into their fit and 
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Mazets etal . ^ ' did not. This leads to somewhat higher radiation 
temperatures than otherwise. Of the published spectra of Mazets and 
coworkers, nearly all can be fit with a spectral form 
dN/dE a ^ exp (-E/E ). The observed range of radiation temperatures 
Eq is quite broad — from ~ 30 keV to more than 1 MeV. Tnat a single 
spectral form provides an excellent fit to the great majority of 
events over such a broad range of radiation temperature is a remarkable 
and significant feature of GRB's and one that GRB models must now 
incorporate. In Figure 4 the distribution of the number of events 
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versus maximum radiation temperature is plotted^ The distribution 
is broad with a maximum at 200-300 keV. 

Recently the first results from the y-ray spectrometer or. the 
Solar Maximum Mission (SMM) have been reported by Share et al . 

In slightly more than a year the instrument has detected 17 
events in the energ\' range 300 keV to 8 MeV. The large vol’ume of 
Nal detectors make this the most sensitive instrument for GRB 
spectroscopy thus far flown. Unf ortunatelv however, the low enercrs’ 
threshold is at ~ 300 keV. Tney find in at least one case (1980 
April 19) that the spectrum is consistent with an E”^’^ power law in 
the .3-8 .MeV interval. Since their energy window is different from 
that of Mazets etal. it is not clear that this conflicts with 
the Mazets observations. It may be, however, that they have found 
evidence for a high energy tail that departs from the thermal 
bremsstrahlung-like spectrum. 
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and 1979 January 13 do appear to display a correlation between 
radiation temperature and intensity* 'Wie former shows a radiation 
temperature thfit peaks well into the event, and the latter shows a 
definite correlation of the temperature with subsidiary peaks in the 
count rate. It is evident that there is a rich variety of behavior 
present in the continuum spectra of GRB's. Some of the implications 
of this behavior will be discussed later in this paper. 

EMISSION AND ABSORPTION LINES IN GAMMA-RAY BURST SPECTRR 

a) Cyclotron Emission and Absorption Features 

During the past year it has become apparent that emission and 
absorption features are common occurrences in the spectra of 
GRB' s ^ ^ The most prevalent type of spectral feature 
that is seen is either an absorption line or broad absorption band 
in the 30-70 keV range^®'^^. Of the ~ 150 events recorded by the 
KONUS experiment, such features are seen in at least 30^*^. The 
proportion of events containing lines may, in fact, be even larger 
since the poor statistics of the smaller events may be masking the 
presence of these features. Two exabmples of "cyclotron" absorption 
features are given in Figure 7a, b^*^. The spectrum of the 1*!79 March 
7 event (Figure 7a) clearly exhibits an absorption feature at ~ 45 
keV having a FWUM of 15-20 keV and an equivalent width of 7 ± 0.6 
keV^*^. The 1978 November 15 event^^ (Figure 7b), on the other hand 
displays a broad absorption laand extending from ~ 100 keV down to at 
least the 30 keV threshold of the instrumnent. Note that there is 
excellent agreement between the spectra measured on the two 
different spacecraft. Evolution of the "cyclotron" absorption 
feature during an event is shown in Figure 8^®* The first 8 sec. of 
the event plotted in the left panel of the figure clearly shows an 
absorption line at 65 keV superimposed on a thermal-bremsstrahlung- 
like continuum spectrum having a radiation temperature of 240 keV. 
The right hand panel shows the subsequent 24 sec where the line 
feature has completely disappeared. Furthermore, the radiation 
temperature has increased to 480 keV. In nearly all of the events 
studied by Mazets et al . ^ ^ the absorption features at^ present 
for less than the total duration of the event. Furthermore, they 
are always strongest during the initial phase and tend to decay away 
more rapidly than the overall burst intensity. In general the 
behavior of the absorption feature is not well correlated with the 
behavior of the continuum spectrum of the burst^^. 

At the present time, only one possible confirmation of the 
presence of cyclotron absorption features in GRB's is known to us. 
Tnis is the event of 1980 April 19 as recorded by the Hard X-Ray 
Burst Spectrometer on S.MM^®. The general character of the 
absorption feature as well as its time variability appear to be 
consistent with the behavior of the events reported by Mazets 
^.10,11^ Tne issue, however, is somewhat confused by the fact that 
a solar origin cannot be ruled out. It should be pointed out that 
solar origin any of the Idazets events is ruled out since directions 
are determined for each event by the KONUS experiment. 








Figure 8 . Spectrum from the 1979 Nov 1 GR3 from the KONUS 
experiment. Left hand panel is the first 3 sec. of the event. 
Right hand panel is the subsequent 24 sec. 


k very striding exaunple of spectral evolution during a burst is 
aiven in Figure 9 where the usual ''cyclotron" absorption feature is 
seen during the early phase of the event superimposed on a 400 keV 
continuum spectrum'*’^. Note, however, that the lowest energ\- point 
at ~ 20 keV departs significantly from the thermal-bremsstrahlung 
representation, indicative of the possible presence of a separate 
soft component. k second spectrum is plotted that was taken during 
a minimum intensity period of the event occuring at ar. intermediate 
time within the burst. The two lowest energ^j’ data points have not 
changed. Tne hard 400 keV spectrum previously present above ~ 70 
keV has, however, now completely disappeared. Only a very soft 
spectram, wh.ich appears as a continuous extension of the unchanged 
low energ\' data points, now remains. Tnis event seem? to have two 
components, a soft one that remains more or less constant and a hard 
one that is time-variable. We shall return to t'^is point later as a 
number of ot'var observations suagest that a two-comoonent monel is 
anpropr iate . 

Tne detailed characteristics of the "cyclotron" features appear 
to hear little or no relationship to the properties of the continuum 
spectram. One can use the Mazets et al . oublished data to 
calculate an effective "FWHK" by dividing the equivalent width by 
the deptn o*' the absorption feature. For the broad band features 
this "FWHM" should be interpreted as a lower limit on the true width 
since the features generally extend below the instrumental threshold. 
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Figure 9 . Spectra 
from the 1979 Mar 
29 GRB from the 
KONvJS experiment. 
Solid line is from 
the initial 
stage. Dashed line 
is from a deep 
minimum between 
burst peaVs. 
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This effective "FWHM" is plotted in Figure 10 versus the continuum 
ra'iation temperature for the Nazets events where data are 
available There is obviously no evidence here for any 
correlation between the line width and the continuum temperature. 
This can be taken as further evidence that two different processes 
or regions are responsible for the line and continuum behavior. 

One final example of a "cyclotron" feature is given in Figure 
This is the 1979 May 26 event which is the only published 



Figure 11 . 
Spectrum of the 
1979 May 26 GRB 
from the KONUS 


experiment 
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example of a "cyclotron" emission feature. It appears at an energy 
of 45 ± 5 keV and has a width of ~ 20 keV. 

The persistent presence of these absorption and emission 
features in the 30-70 keV range of GRB spectra is strong evidence 
for a neutron star origin. In the intense (10^^-10^^ gauss) 
magnetic fields that are expected to be present near the surfaces of 
neutron stars, the cyclotron frequency corresponds to photons in the 
hard X-ray range. The fact that such lines exist superimposed on a 
very hot continuum spectrum poses a problem. To produce a cyclotron 
absorption line requires that the electron temperature be smaller 
than the energy’ of the line (first Landau level). If the tempera- 
ture were comparable to or higher than this, one would expect to see 
the line in emission. The one example of "cyclotron" emission (197° 
Kay 26) may, in fact, be a case where the temperature in the 
emission region is higher than the first Landau level. The width of 
the line (in some cases ~ 20 keV) imposes a further constraint on 
the temperature of the absorbing region. The Doppler broadening of 
the line is given by Trumper^*^: 
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~ ■ /8tn2 — j cos 6 
me 

where 9 Is the angle between the line of sight to the observer and 
the magnetic field. Taking ~ ■ 0.3 and 0 » 45° gives kT ■ 16 keV 
for the electron temperature. To account for these observations it 
appears that one must postulate that the cyclotron absorption takes 
place in a relatively cool (< 20 keV) layer that overlies the hot 
region which produces the continuum spectrum. 

We return now to the question raised in the previous chapter of 
whether optically thin thermal bremsstrahlung can be a realistic 
physical description of the radiation process during a GRB. The 
presence of cyclotron features in the 30-70keV range implies that 
10-10 gauss fields are present in the emission region. It is 
then extremely unlikely that free-free transitions will be the 
dominant mode of radiative transport. It is much more likely in 
this intense field that synchrotron radiation and related processes 
will dominate. 

b) Enission Lines > 400 keV 

Bnission lines > 400 keV in GRB spectra have been reported by 
Mazets et al .^^*'*’’ and Teegarden and Cline"* The former, as 
mentioned earlier used Nal scintillators, whereas the latter 
employed a radiatively- cooled Germanium detector. This instrument, 
launched on the lSEE-3 spacecraft, had a resolution of 8 keV at the 
time of launch. Unfortunately, after only 4 months of operation , 
an electronics failure disabled that portion of the memory in which 
GRB spectra were stored. As a result, high-resolution Ge spectra 
are available for only two events. 

An example of an emission line at ~ 400 keV is given in Figure 
12 for the event or 1978 September 18*®. The line appears quite 
pronounced and is superimposed on a 185 keV thermal-bremsstrahlung- 
like continuum. As was the case with the **cy''lotron" lines, the 400 
keV features do not appear to bear any correlation with the 
continuum radiation temperature. Of the 150 events recorded by the 
XONUS experiment, 11 display emission features in the 400-460 keV 
range. These featur have been interpreted as redshifted 
annihilation radiation. Redshifts in the 15-20 percent range are 
expected for radiation produced near the surface of a 1 neutron 

star. This would shift the 511 keV annihilation line into the 400- 
460 keV range. There are, unfortunately, a number of factors that 
complicate this interpretation. First, for high electron 
temperatures the energy of the line will be blueshifted since all of 
the energy of the electrons (kinetic + rest mass) is converted into 
photon energy^® It is not clear, however, that electron 
temperatures high enough for this effect to be significant are 
consistent with the observations. In addition, the possibility of 
stimulated e*-e~ annihilation exists if one is dealing with a cool 
dense electron plasma (Ramaty, private communication). It turns out 
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Figure 12 . Spectrum of the 
1978 Sept 18 GRB from the 
KONUS experiment.^® 



Figure 13 . Spectrum 
of the 1979 June 22 
GRB from the KONUS 
experiment^® 


that the energy dependence of the photon absorption coefficient is 
such that a peak in the 400-460 keV range can be produced. 

Figure 13 is an example of an event where both a “cyclotron" 
absorption feature and a 460 keV emission line are observed^® 
simultaneously. There is only one other such example that has been 
published by Mazets and coworkers This is the event of 1979 
November 1 which is reproduced in Figure 14. This single figure 
displays much of the rich phenomenology tnat has manifested itself 
in the spectra of GRB's. The left panel of the figure shows the 
^irst 8 sec of the event where both a 65 keV absorption line and a 
bread emission feature in the 350-650 keV interval are seen. These 
lines lie on top of a 280 keV continuum. A profile of this emission 
hand with the continuum subtracted is also shown. T\’pically, the ~ 
400 keV emission features have widths of ~ 250 keV, in the right 
hand panel two spectra in the later portion { resnectively ~ 16 sec 
and ~ 24 sec after the onset) of the event are displayed. Both the 
absorption and emission lines have disappeared. In fact, the 
typical behavior for the ~ 400 keV features is appearance at the 
onset of the event and a lifetime shorter than that for the total 
intensity. Such behavior, as was seen earlier, is also commonplace 
’•ith the "cyclotron" features. During the intermediate portion of 
the event, the radiation temperature rises to a value of 800 keV and 
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Figure 14 . Spictra of the 1979 Nov 1 GRB from the KONUS 
experiment’' "I . (1-2) first 8 sec of the event; (4) 4th 4-sec 

interval; (6-8) 6^^, 7-^ and 8^^ 4-sec intervals. 


then decays to a value of 440 )ceV in the latter portion. Curiously 
enough, the disappearance of the lines is accompanied by a heating 
of the continuum temperature. The general, lack of correlation 
between the " 400 keV emission features and the continuum 
temperature can be taken as further evidence for the existence of 
two components in the GRB source. Furthermore, the observed 
broadening (~ 250 keV) of the annihilation line requires that the 
electron temperature be low. If the line width were entirely due to 
Doppler broadening, we would have kT 2 15 keV. Furthermore, as 
pointed out by Daugherty and Bussard^^ and Mazets et al . ^ ^ > the 
radiation produced by pair annihilation in strong magnetic fields 
will be broadened by the presence of the field itself. In fact, 
this broadening mechanism in a field of 5 x 10^^ gauss could account 
for the full observed width of the 400 keV lines. 

Ramaty and coworkers^^ ' have proposed a model for the 1979 
Karch 5 event where the radiation originates in a thin (~ 10”^ cm) 
surface layer that overlies a hot MeV plasma. Tney point out that 
electron cooling by syrchroton radiation in a > gauss field 

will generally take place more rapidly than e e“ annihilation. In 
their model the continuum spectrum is therefore produced by the 
s^mchroton emission from the rapidly cooling electrons and positrons 
which subsequently annihilate and produce observable lines whose 
width is characteristic of the low temperature in this thin surface 
layer. The 1979 March 5 event is unique in a number of ways: (1) 

its luminosity is brighter than any other recorded event, (2) eight- 
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second oscillations were observed superimposed on a slowly decaying 
tail, (3) it has been identified with a supernova remnant, N49, in 
the Large Magellanic Cloud and (4) its continuum spectrum is very 
soft. Ticse unique features may point to different physical 
conditions at the source of this event. It is nonetheless true that 
the Ramaty et al .^^*^^ model for the 1979 March 5 event may be 
applicable to other GRB's. 

A cool layer that is optically thin to Compton scattering 
appears to be required to explain the spectral behavior that has 
been recently observed. The fact that in many events ~ 1 MeV 

photons are present would imply that Y - Y pair production in the i 

region immediately above the absorption layer is not an important 

process. Schmidt^^ has used this argument to place a limit on the 

distances of GRB sources. He assumed an emitting region of 10- cm 

extent and obtained a distance limit of ^ 2 kpc for events in which 

MeV photons are present. The weight of recent evidence is towards a 

much smaller size (~ 10^ cm) for the emission region which in turn 

may lead to a smaller distance limit for the burst sources. The 

problem is, however, congslicated by the strong energy and angular 

dependence of the cross-sections^^ eo tha^ detailed calculations 

will be necessary to establish a true distance limit. 

Figuve 15 shows the spectrum of the 1978 November 19 GRB as 
measured by the Germanium spectrometer onboard the ISEE-3 
spacecraft^ A broad marginally significant feature is present 
at 400 keV that is consistent with the spectral feature reported 
by riazets et al .^^ during this same event. In addition, the ISEE-3 
data show a second line at an energy of 740 keV having a width of ^ 

40 keV. This line is signif .cant at the 99 percent confidence 
level^^. Teegarden and Cl'-'o^^ pointed out that the energy of this 
line is consistent with that expected if the 847 keV first excited 
level of iron is redshifted by the same amount as the 400 keV 
feature. It is perhaps significant that this event has the hardest 
spectrvsn of any thus far recorded and is one of the few that 
apparently cannot be fit with a thermal-bremsstrahlung function. 

The continuum spactrum for energies ~ 200 keV-2 MeV best fit with an 
power law^^. 

If the assumption of a redshifted iron line is correct- then 
conditions at the source are probably such that other nuclear levels 
can be excited. ’<»e do not a priori know the con^osition of the 
source region. If it is pure iron, then we would see only higher 
levels of Fe nuclei., but if it ia, for example, the sane as the 
'olar system composition, then excited levels of other nuclei would 
probably produce significant radiation. Ramaty et al -^^ have 
calculated the production of gamma-rays by cosm;.c rays interacting 
with the Interstellar medium. Vte have normalized their iron line 
intensity to that of our 740 keV line and then plotted the most 
prominent lines > 740 keV normalized by the same factor and 
redshifted by the same amount as the Fe line. The Ramaty etal. 
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yjqure 15 . Spectrum of the 1978 Nov 19 event from the lSEE-3 Gamma 
Ray Bunt Spectrometer^ 

calculations^^ that we used assumed an E"^ power law proton spectrum 
which may not be a good representation of ti>e spectrum at the site 
of the GRB, The relative line strengths, however, r.re not strongly 
dependent on the shape of the proton spectrvm but rather on the 
relative abundances and excitation cross-sections. It is evident in 
figure 15 that these higher energy lines could make a sign^icant 
contribution to the > 1 MeV portion of the spectrum and that they 
might, in fart, account for the unusual hardness o* the spsetrum of 
this event. 

Katz (1931) has put forth an alternate explanation for the ’’40 
keV line. He has pointed out that for sufficiently high magnetic 
fields (~ lO^^qausr) single photon annihilation will become the 
dominant mode. In this mode a single photon with energy 1.022 MeV 
is emitted so that the ratio of the two lines should in this sinple- 
mlnded picture be a factor of twr The actual ratio in the data are 
- 1.76, but the peaks are broad e..ough so that a factor of two 
cannot be ruled out. Ramaty et al . ^ ' have pointed out that a 
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mechanism exi'ts which can effectively shift the annihilation peaks 
toward lower e lergies and in general the one- and two-photon peaks 
will not be shifted by the same amount. It therefore remains 
unclear whether this can serve as an adequate explanation of the 
spectral features in the 1978 November 19 event. 

SUMMARY AND CONCLUSIONS 

The most important characteristics of GRB spectral behavior are 
summarized as follows: 

1 ) The continuum spectrum is well represented by a thermal 

brerasstrahlung function '■'ith kT ranging between 30 keV and — 1 

MeV. Possible evidence xists for a power-law high energy 
tail. 

2) Significant evolution of the continuum sprectrum during the 
event is commonplace. The most typical behavior is a higher 
temperature during the early part followed by "cooling" or 
softening of the spectrum. There are, however, many examples 
of events where the temperature appears to be correlated with 
the intensity which may in turn possess complex structure. 

3) "Cyclotvon" absorption features in the 30-70 keV range are 
frequent, if not typical, occurrences in GRB's. At least 30 
out of the 150 events recorded by the KONUS experiment display 
such features. Line widths vary from as small as ~ 20 keV (the 
probable instrumental limit) to broad absorption bands 
extending over more than 80 keV. "Cyclotron" features are 
normally present during the early part of the event and 
disappear before the event is finished. Little or no 
correlation exists between the continuum radiation temperature 
and the properties of the absorption lines. 

4) Emission lines in the energy interval 400-460 keV are present 

in at least 7% of the GRB's. The features apparently are broad 
(AE ~ 200 keV) and uncorrelated with the continuum 

temperature. As with the cyclotron features they appear at the 
onset and disappear before the event is over. 

5) A narrow (" 40keV) line has beer, reported in the 19 78 November 

19 event at 740 keV. This event has the hardest spectrum of 

any thus far recorded and is one of the few that apparently 

cannot be fit by the thermal-brehmsstrahlung expression. 

Tne principal conclusions that can be drawn from these observations 
are summarized as follows: 

1) The "cyclotron" and 400-460 features car be ,aken as strong 

evidence for a neutron star origin for GRB's. Tne 400-450 keV 
lines are interpreted as redshifted annihilation radiation. 

The redshift is consistent with that expected from a 1-1.5 M*, 
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neutron star. 

2) The presence of absorption features <70 keV and ~ 400 keV lines 
with AE/g ~ 0*5 appears to require the existence of a cool 
layer overlying the hot plasma in which the continuum is 
produced. Rapid synchrotron cooling is a possible mechanism 
for the production of such a layer. 

3 ) The presence of ~ 1 MeV photons in many events requires that 
the distance to the source D ^ 1 kpc. At larger distances y-y 
absorption would strongly suppress the high energy photons. 

4) Optically thin thermal-brehmsstrahlung is probably not a viable 
description of the radiation mechanism at the source. In the 
strong magnetic fields in the source region synchrotron 
emission and absorption are almost certainly much stronger 
processes. Any theory, however, must account for the fact that 
the continuum spectra of nearly all GRB's are extremely well 
fit by a thermal-brehmsstrahlung expression over a wide range 
(30 keV - 2 MeV) in energy. 

5) The 740 keV line could be explained either by a redshifted iron 
line or single-photon annihilation. 
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ABSTRACT 

The ISEE-3 ganma-ray spectroneter has provided nearly 
continuous monitoring of the cosmic gamma-ray background in the 
energy regime 125 keV to 6.5 MeV since launch of the satellite In 
August 1978. The data has been analyzed for possible variations in 
the cosmic 511 keV line flux. The detector Is an unshielded, 
radlat ively-cooled, hlgh-purlty germanium crystal with a sensitive 
volume 33 cm'*. At energies > 100 keV the detector has a fleld-of- 
vlew of > 2n steradlans which Includes the galactic center. The 
detector resolution v;as degraded continuously by exposure to high- 
energy charged particles In the Interplanetary space environment; 
consequently, only the first 500 day sample of data Is usable. 
Proton and electron flux rates available from a cosmic-ray 
experiment aboard ISEK-3 were used to eliminate all periods during 
which solar energetic particles contributed significantly to our 
counting rate. Ue are able to place upper limits on the 
variability of th(^ c<^Y>lc 511 keV flux of < 2.2 x 10~ and < 1.3 x 
10 photons cm *■ s for = 20 day and 60 day integrations, 
respectively, for the period October 1978 to February 1980. 
Comparison is made with the HEAO-3 511 keV observations of the 
galactic center. 

INTRODUCTION 

During the past decade several observations of the galactli 
center region have produced positive detectlcws of emission in the 
vicinity of the 511 keV annihilation, line ** ’ ’ ’ Th/f 

observed fluxes range from 0.8 x lO”"* up to 4.2 x 10~^ photons cm”*" 
s . The first observations from a satellite-borne experiment 
appear to support the Interpretation that the emission may 
originate from a time-varying ‘point*’ source. Since the available 
data have been sporadic, It would be desirable to have a continuous 
record of observations over an extended period. We report here 
upper limits to the variation of galactic 511 keV flux as observed 
by the ISEE-3 gamma-ray spectrometer for 500 days of nearly 
continuous observation. 


INSTRUMENTATION 

The ISEE-3 (Third International Sun-Earth Explorer) spacecraft 
Is In a halo orbit about the Lagranglan point on the earth-sun line 


230 earth radii (4.5 light-seconds) inward towards the sun. This 
location Is particularly advantageous; the time-varying exposure to 
trapped radiation within the magnetosphere and the earth-albedo 
problem are completely avoided. 

The spectrometer is approximately aligned with the spacecraft 
spin axis, that is, normal to the ecliptic plane. The Instrument 
is recessed within the lower body of the spacecraft, preventing 
direct sunlight from destroying the thermal performance of two 
nested radlatl vely-cooled stages. The detector is housed In the 
inner stage. The outer stage Is conically shaped to define a I2b° 
(thernsjl) f leld-of-vlew for the Inner stage. This f leld-of-view 
does not contain the sun, earth, or any significant portion of the 
spacecraft . 

The detector is an unshielded p-type co-axlal hlgh-purity 
germanium crystal, 4 cm dlam. x 3 cm depth, with an active volume 
of 33 cm’’. Pre-launch tests measured an energy resolution of 3-3.5 
keV at 1 MeV. A higher temperature (130 K) than the predicted 
equilibrium temperature of 100 K for the crystal was reached when 
the spacecraft was placed In orbit. Consequently, the Inj^tlal 
operating resolution of a calibration line at 570 keV from a *■ 
source was ~ 8 keV (FWHM). 

OBSERVATIONS AND ANALYSIS 

The experiment was Intended to monitor the celestial sphere 
for gnmma-ray bursts, providing information on any narrow spectral 
features within the energy range 125 keV to 6.5 MeV. In the 

absence of a burst trigger, a b-ackground mode operates, permitting 
spectral analysis of the diffuse gamma-ray background radiation. 
The background mode has functioned successfully since September 
1978. Because of the extremely low bit rate allocated to this 

experiment, the Instrument analyses only one photon of ' 550 

detected In the background mode. Hence, in order to accumvilate a 
sufficient number of counts to constitute an observation, the 
spectra must be accumulated for a minimum of several days. 

F.nergy calibration wa« achieved by monitoring t» o background 
lines at 570 keV and 1064 keV. For the purpose of rebinnlng counts 
into fiducial channels, it was assumed that the gain was a linear 
function of energy. The positions and widths of the lines lu 
accumulated spectra were fit both manually and by an automatic 
gaussian line fitting routine. The gain (channels/keV) decreased 
very nearly linearly by about lOX during the first 650 days of 

experiment operation. Also, significant deterioration of ti>e 

energy resolution occurred: ti>e 570 keV line broadt^ned fr>m the 

initial 8 keV FWHM to ' 50 keV in 600 days (Figure 1). After this 
period, tlie 511 keV and 570 keV lines began to overlap 
8 1 gnl t leant ly . The relative magnitudes of the gain and resolution 
changes .i.cc consistent with a decreasing charge collection 
efficiency, evidently caused by detet-tar radiation damage from 
li 1 gh -c ne r gy .'barged particles and fast neutrons. 

Time variation In the Interplanetary charged particle flux due 
to solar modulation and other effects will produce variations in 



Fig. 1.- Resolution of 570 keV calibration line for 600 days from 
start of experiment. 

the detector background that must be evaluated. Particle fluxes 
for several different proton energy ranges and for electrons were 
Independently recorded by a cosmic-ray experlme|t aboard ISEE-3 and 
there data were kindly made available to us . Those particles 
(energies = 14 MeV to 18 MeV) which deposit energy in the detector 
within the PHA energy range originate primarily from solar 
flares. Therefore, to calculate the actual particle rate 
contribution, a canonical solar flare energetic particle spectrum 
proportional to E * was assumed and the required rates 
extrapolated from the available data. Periods for which the 
calculated particle rate exceeded one-tenth of the random error for 
given spectral accumulation Intervals were eliminated from the data 
base. From late 1978 to 1980 the sun was especially active. 
Consequently, the requirement for non-contamlnated spectra resulted 
in eliminating up to 60% of the data. 

Extraction of the astrophyslcal 511 keV flux from the 
“cleaned" data was complicated by two effects. First, a large 
fraction of the detected flux Is instrumental in origin and Its 
contribution is difficult to estimate. Furthermore, the capability 
to point the spectrometer does not exist; hence an on-source minus 
off-source type of observation cannot be conducted. At energies 
> 100 keV, the cooler stages and the spacecraft well into which the 
spectrometer Is recessed are virtually transparent, affording the 
detector an effective f leld-of-vlew > 2v steradians. The galactic 
center region, from which a variable 511 keV flux may emanate, Is 
constantly within the detector f leld-of-vlew. Assuming the 
galactic center is the most intense 511 keV source, for a galactic 
center observation It is sufficient to distinguish a variation in 
the total 511 keV flux. This subject will be addressed In the 
Discussion section of this paper. A second complication results 
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from the broadening Instrumental resolution due to detector 
exposure to particle radiation. Because of the low statistical 
significance of the accumulated spectra, it was not acceptable to 
employ a straightforward technique of fitting a gaussian line to 
the 511 keV feature and subtracting a continuum. Instead, to 
accumulate counts of = 511 keV energy, an automatic "window" 
spectra program was employed which appropriately gain-adjusted the 
spectra and summed counts in four energy windows of widths 38, 50, 
60, and 70 keV, centered at 511 keV. This treatment results in a 
systematic decrease of counts with time in a given energy window 
due to the broadening resolution. Therefore, a reduction technique 
was employed for each accumulated 511 keV window spectrum in 
which 1) an adjusted continuum computed from a window spectrum of 
width 40 keV, centered at 450 keV, was subtracted from the line 
plus continuum spectrum (the continuum above 511 keV was not used 
due to the presense of the strong calibration line at 570 keV); 2) 
a gaussian— shaped line which broadened with time commensurate with 
the observed broadening of the calibration lines was assumed for 
the 511 keV line, and normalized model curves were constructed 
which accounted for the systematic loss of counts from the energy 
windows; 3) the observed count rates and the models were 
subtracted to yield relative 511 keV count rate variations as a 
function of time. The count rates were then converted to fluxes by 
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Fig. 2 .- One sigma upper limits on variation of 511 keV flux. The 
Integration Intervals (live times) are approximately 20 days (top) 
and 60 days (bottom). The data are plotted with respect to the mean 
values (dotted lines). 
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dividing by the effective area of the detector. The wider energy 
windows have an Inherently poorer sensitivity due to the inclusion 
of aore background counts. For the narrowest (38 keV) window, the 
sensitivity decreases more rapidly with time as the Sll keV line 
progressively broadens, leaving less counts within the static 
window. The sensitivity Is effected more drastically by the use of 
wider energy windows. The resulting fluxes for the 38 keV window 
are shown In Figure 2 for a 20 day and ~ 60 day integrations (live 
times!. The two sigma upper limits are 2.2 x 10“^ (20 day) and 1.3 
X 10 (60 day) photons cm s on the 511 keV flux variation from 
> 2 ti steradlans for the period September 1978 to February 1980. 
During the same period, the (> 60 MeV) particle flux determined 
from the ISEE*-3 cosmic-ray experiment was practically constant 
(± < 4Z) for the intervals of retained data. The instrumental 
background results primarily from Illumination of the spacecraft 
and activation of the detector crystal by this radiation. The 
percentage particle flux variation for any given spectral 
accumulation Interval accounts for a maximum of less than one- 
fourth of the quoted limits on the 511 keV flux variation. No 
correlation is apparent between the (> 60 MeV) particle rate and 
our final results. 


DISCUSSION 

The galactic center observations performed with the HEAO-3 
gamm.i-ray spectrometer support the interpretation thac a time- 
varying source of 511 keV emission lies pesc the galactic center 
Several balloon experiments ’ ’ ’ ’ • Indicate that the 

source flux may range from 0.8 x 10 up to 4.2 x 10 photons cm”^ 
8~ . The lSEE-3 upper limits for the 511 keV flux may be examined 
in two ways. First, since the spectrometer f leld-of-vlew includes 
more than half of the galactic plane, the ^ 20 day observations 
exclude transient galactic sources (i > one month) with fluxes > 
2.2 X 10 photons cm s (952 confidence). If sources similar to 
the galactic center source were distributed uniformly thoughout the 
galaxy, a source within 2-3 kpc could have been detected by the 
ISEE-3 spectrometer. The lack of detection Is evidence that the 
galactic center source is unique within our galaxy. Second, the - 
125 day observations (60 day live^time) constrain ^he peak galactic 
center source flux to < 1.3 x 10 photons cm s during a period 
of 500 days from October 1978 to February 1980. 

Wc thank T. von Rosenvinge for providing us with particle data 
from an ISEE-3 cosmic-ray experiment. We also thank W. Paclesas 
for the automatic gausslan line fitting routine and for helpful 
discussions. J. P. Norris acknowledges support from a NASA Student 
Researchers Grant. 

From a dissertation to be submitted to the Graduate School, 
University of Maryland, by Jay P. Norris, In partial fulfillment of 
the requirements for the Ph. D. degree In astronomy. 
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ABSTRACT 

A short but intense y -ray burst was observed on 1980 April 19 
at 01^ '9° 46® with the Hard X-ray Burst Spectrometer on the 

Solar Maximum Mission. The event lasted for a total of 10 s and 
consisted basically of three 200 - 300 ms wide spikes with signif- 
icant variability on a time scale of 20 - 40 ms. The photon number 
spectrum integrated over the impulsive part of the event fits a 
thermal brensstrahlung function with a temperature of (330 ± 70) keV 
at energies from 151 to 487 keV. At lower energies the data points 
lie considerably below this function suggesting a broad absorption 
feature extending down to < 28 keV, the lowest energy measured. The 
upper energy of this absorption feature varies between 100 and 150 
keV on a time scale of <0.5 s. This event is interpreted as a 
typical gamma ray burst although it is still remotely possible that 
it is of solar origin. The spectral features and their variability 
are interpreted in terms of electron interactions at the cyclotron 
resonance frequency in magnetic fields of 10^^ - 10^^ gauss close 
to the surface of a neutron star. 

INTRODUCTION 


Absorption and emission features in the high energy X-ray 
spectra of compact . astrophysical objects are believed to result 
from magnetic bremsstrahlung in intense fields of the order of 
10^^ - 10^^ gauss. At these intense fie.ld strengths the plasma cy- 
clotron frequencies are in the 10 to 100 keV energy range. Emiss- 
ion and absorption at such frequencies create spectral features 
that can be used to determine magnetic and plasma conditions in 
the vicinity of the compact object. 

The first observation of such a feature was made by Trunper 
et al . ^ in the spectrum of Her X-1. This feature, now confirmed by 
Gruber et al .^, can be interpreted either as a cyclotron emission 
line at - 58 keV or as an absorption dip at - 42 keV. Such emiss- 
ion and absorption features have also been detected in the spectra 
of some gamma-ray bursts by Mazets etal. They observed absorpt- 
ion features extending below - 100 keV in the spectra of 20 out of 
- 150 gamma ray bursts de*-ected on the Venera 11 and 12 space probes. 
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and they observed an emission feature at 45 keV in the spectrum of 
one burst. They obtained spectra every 4s throughout the events and 
in sone bursts observed large variations in the observed features. 
Thtir observations of absorption and emission lines interpreted as 
cyclotron features provide the most convincing evidence to date that 
gamma ray bursts are produced by neutron stars. 

In this letter we present high energy X-ray spectral data for a 
gamma-ray burst with over an order of magnitude finer time resolution 
than that obtained by Mazets et al .^ The measured photon number 
spectrum can be interpreted as showing a broad absorption feature at 
energies below 100 keV simi'.ar to the features reported by Mazets 
c_c for other events. Tiiese obser'’ations support the idea 

that gamma-ray bursts are produced in the vicinity of neutron stars. 

OBSERVATIONS 

The gamma ray burst was detected on I960 April 19 at 01^ 19*“ 46® 
with the Hard X-ray Burst Spectrometer (HXRBS) on the Solar Maximum 
Mission (SMM). The HXRBS is an actively shielded and collimated 
CsI(Na) scintillation spectrometer designed primarily to detect high 
energy X-rays from solar flares^. Although the instrument is point- 
ed at the sun at all times, any off-axis source is detectable within 
a 40° FWHM circular field of view. Outside this field of view the 
central detector is shielded by a minimum of 1" of CsI(Na) operated in 
anti-coincidence. 

The arrival direction of the burst was determined to be '> 18" 
from the sun using the detection times at SMM, Vela 5 and Vela 6 in 
earth orbit, at ISEE-3 located 4.5 light seconds away in the solar 
direction, and at the Pioneer-Venus Orbiter^. There is still a 
remote possibility, however, that the source was of solar or near- 
solar origin. Nevertheless, it is clear that the source was well 
within the field of view of the HXRBS, and no significant spectral 
distortions should have occurred. We also know that, at the time of 
the event, the Earth was well out of the HXRBS field of view, the 
horizon being > 50“ off axis. Thus, fluorescence from the earth's 
atmosphere could not have contributed to the source spectrum measured 
during the burst. 

The time profiles of the event in several different energy bins 
covered by the HXRBS are shown in Figures I and 2. The event is 
typical of many gamma-ray bursts and consists basically of three 
sharp, 200 - 300 as wide spikes separated by 500 - 600 ms. Signlf- 
i<'ant variations in flux are evident o* time scales as short as 20 - 
40 ms. A gradual flux Increase precedes the first spike by 1 to 2 s 
and a general decay is evident for several seconds after the third 
spike. The fluctuations in the decay phase of the event do not show 
significant power at any frequency from 0.5 to 10 Hz. 

The time profiles in the different energy ranges plotted in 
Figure 2 show clearly that the spectrum changes significantly through- 
out the event. The first spike is barely visible above a general 
plat-nu at energies below 126 keV but it shows up clearly in the 126 
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FIGURE 1 : Tine profiles of the ganna-ray burst in different energy 

ranges and with different tine rer.olutions. (a) Counting rate in 
the central detector for the energy range fron 28 to 487 keV with a 
tine resolution of 128 ms. The tine intervals used for the energy 
spectra in Figures 3, 4 and 3 are indicated. (b) Counting rate in 
the central detector for the same energy range but with a tine res- 
olution of 20 ms. This data has not been corrected for instrunent 
dead time, which changed from 5X before and after the burst to 
'' 157. at the tine of the biggest spike in Interval ^<4. (c) Count- 

ing rate in the CsI(Na) shield crystal for all energies 2QU keV 
with a time resolution of 64 ns. The spike at 01^ I9n 53s and all 
fluctuations below •- 1200 counts are variations in the back- 

ground rate. 


to 151 keV and the 203 to 232 keV energy ranges and also in the 
shield data at ^200 keV shown in Figure 1. The second spike is 
clearly resolved only in the shield date. The third spike is also 
not well resolved at energies below 78 keV but shows up strongly 
at higher energies. 

The photon number spectrum for the impulsive part of the 
event (intervals 2, 3 and 4 as Indii-ated In Figure 1) is plotted 
in Figure 3. The spectral points were obtained iteratively by de- 
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FIGURE 2 : Tine profiles of Che 
ganna-ray burst for the indi- 
cated energy ranges all with a 
cine resolution of 128 ms. 


convolving the observed 15-chann- 
el counting race data through Che 
instrunenc response natrix for an 
assumed incident spectrum equal 
to Che least-squared fit exponen- 
tial function indicated in Figure 
3. The instrunenc response matrix 
was computed for different inci- 
dent spectral shapes using the 
energy calibration determined for 
pre-launch data and in-orbit data 
from the the on-board ' 
source^. The computations take 
into account the detector effi- 
ciency and resolution, I-cscape 
and Compton scattering using the 
results of pre-launch ‘..>urce cal- 
ibrations and Monte Carlo simu- 
lations. To determine the accu- 
racy of the calculations, compar- 
isons were made of HXRBS spectra 
for several solar flares with 
spectra of the same flares ob- 
tained with the Gamma Ray Experi- 
ment (GRE) on SJIM (Ryan, private 
communication), and with the hard 
X-ray detector on ISEE-3 (Kane, 
private connunication) , These 
comparisons suggest that syste- 
matic uncertainties are less than 
± 30% for any of the 15 channels 
although all of the flare spectra 
are considerably steeper (E"'^) 
than the spectrum of the gamma ray 
burst. An additional systematic 
uncertainty arises from the prob- 
able 18* off axis position of the 
source. This results in a reduc- 
tion in sensitive area ', / a factor 
of 2 below the 68.6 cm- for an 
on-axis source, with the conse- 
quent increase in the calculated 
photon fluxes by the same factor. 
As Indicated in Figure 3, an expo- 
nential spectrum of the form 

photon flux - A exp(-E/EQ) 

is an acceptable fit to all the 
data points. Here E is the photon 
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F.TGURE 3 : Photon number spectrum for the intervals 2, 3 srtd A 
indicated in Figure 1 (impulsive part of the event). The error 
bars are ± lo based on the square-root of the observed number 
of counts in each energy interval. The solid curve represents the 
thermal brensstrahlung function with kT ■ 330 keV, which gave the 
least-squares fit to the data points above 151 UeV (channels 7-15). 
The broken curve vcpresents the exponential function giving the 
best fit to ill the data points. 


energy and A and are free parameters detenained in the fitting 
procedure. The value of obtained is 159 ± 8 keV and 1* 

13.5 for 13 degrees of fi-;;ion with the counting rate statistics 
alone used to weight the points. A double power-law function also 
fits the data equally well with a spectral index of -0.5 i 0.1 at 
energies below 300 keV and -A.O i i.O at higher energies. The 
apparent dip in the spectrum in channel 3 is only significant at 
t the 2.2o level below the fitted exponential function. It nay 

result from an error in the value used for the thickness of the 
dead layer on the Csl central crystal^. This dead layer, est- 
i inated to be 5 mils thick from measurements of the ratio of t.ie 22 

I • keV to the 88 keV photon fluxes from results in a dip in 

I the count rate spectrum at energies aoove the K edges of cesium 

I and iodine (33 and 36 keV respectively). Consequently, we would 

expect any inaccuracy in the assumed thickness of this dead layer 
to result in a dip in channel 2 although such a dip could appear 
in channel 3 if the incident spectrum had a steep low energy com- 
ponent. No such, dip in channel 2 or 3 has been observed in any 
solac flare spectrum. 
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Following Mazets _et and Gilnan et al.Q, we have also used 

the function expected for brensstrahlung from an optically thin 
plasma at a temperature, T. This function of the form 

photon flux “ B exp(-E/k.T) 

where k is Boltzman's constant, a results from the temperature-aver- 
aged Gaunt factor, and B and I are free parameters determined in the 
fitting procedure. The value of o changes from 0.48 to 0.148 as kT 
varies from 10 keV to 1 MeV but for comparison with the results 
of Mazets ^ a^. we have taken a = 0 in all of our analysis. 

Using the correct value of a would increase the significance of the 
absorption feature and also increase the least-squares fit temp- 
erature. An acceptable fit of this thermal brensstrahlung function 
to all 15 data points is not possible since it can never be flatter 
than However, an acceptable fit can be obtained to the 

data for channels 7 to 15 corresponding to the energy range from 151 
to 487 keV and this is shown in Figure 3. The best fit value of kT 
is 330 ± 70 keV. At energies below 151 keV the data points lie con- 
siderably belov; the level of this function, at least a factor of 2 
lower at 60 keV, suggesting a broad absorption feature extending 
down to < 28 keV, the lowest energy measured. For the purposes of 
investigating spectral changes with time, the event was divided in- 
to the eight contiguous intervals indicated in Figure 1. The spec- 
tral data shu-.’n in Figure 4 was obtained assuming the best-fit 

exponential spectrum tc rJi the data points in each interval. 
The thermal bremsstrahlurg function was then fitted to as many of 
the higher energy points as possible while maintaining an acceptable 
value of x“* Th^ resulting curves are plotted in Figure 4. 

The spectra for the eight intervals show a general softening 
with tine. The temperature obtained from the fits changes from 

400 keV for interval //I to 70 ± 20 keV for interval #8. The ab- 
sorption feature appears in all intervals except and possibly //5. 
In intervals 1, 2 and 4 the absorption feature extends from at least 

150 keV down to < 28 keV, the lowest energy measured. In intervals 

3, 6 and 7 the absorption appears to extend only to ~ 100 keV or 
less. The apparent emission features at 330 keV in interval #2 and 
260 keV in interval #3 are only significant at the - 2o level and, 
in any case, they are narrower than the instrument could allow. 
The spectrum or interval #7 can be interpreted as containing an 
emission feature at 100 keV with a high temperature continuum. 
The lower temperature fit (kT •• 50 ± 20 keV) with aUzorption below 
- 100 keV as indicated in Figure 4 appears more likely, hovjever, 
in view of the normal "cooling" observed by Mazets ^ many 

gamma ray bursts, and the low temperature continuum (kT ■ 70 i 20 
keV) measured in interval #8. 



FIGURE 4: Same as Figure 3 but for the individual intervals 1-8 

indicated in Figure 1. The curves represent the thermal brems- 
strahlung function giving the least-squares fit to as many of 
the higher energy data points as possible while still obtaining 
an acceptable value of The values of kT in keV for the 

indicatei curves are 406, 229, 510 and 302 for intervals 1, 2, 

3 and 4, and 375, 229, 52, and ^3 for intervals 5, 6, 7, and 8 
respectively. 

DISCUSSION 


It is now almost certain that the April 19 burst was ot non- 
solar origin. The timing information from five separate spacecraft 
gives an error box ^ 18° from the sun. The rapid flux variations, 
the short duration of the event, and its extremely hard spectrum are 
all typical cf gamma ray bursts. Furthermore, the sun was very quiet 
at this time with no reported Ha flares. A small soft X-ray event 
was in progress at the C2 level as observed by the GOES 2 and GOES 
3 satellites (Donnelly, private communication), and a small GRF 
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radio burst was observed at 3750 MHz fron NOAA region #2396 near 
the east limb at the same tine (Enone, private conraunication) , 
The SMM pointed instruments were observing region #2389 near the 
west limb and saw nothing. A coronal transient was observed off 
the N.W, limb by the SMM Coronagraph/Polariueter in this tine frame 
(Sawyer, private communication) but was probably not connected with 
the X-ray event. Thus, it is clear that, if the hard X-ray event 
was of solar or near-solar origin, it would constitute a new type 
of event never previously detected. Consequently, we prefer to 
assume that the event is a ganna ray burst and will interpret the 
data on that basis. 

The photon number spectrun obtained for this event cannot be 
interpreted unambiguously. Observations of gamma rays to energies 
in excess of 5 MeV in this e”ent^*^ show that the simple exponential 
fit to the HXRBS data alone is probably not correct at higher ener- 
gies. A two-power-law fit could, however, represent the data. It 
is unclear at present if the thermal brensstrahlung function could 
also fit all the high energy data above ~ 150 keV. In any case, 
there are many difficulties associated with the assumption that 
the emission is thermal bremsstrahlung from an optically thin hot 
plasma^. Radiation other than from free-free electron transitions 
would predominate fron the expected nonequilibrium, relativistic 
and strongly magnetized plasma near the surface of a neutron star. 
Also, this plasma, while optically thin for free-free transitions, 
would be optically thick for Compton scattering with consequent 
strong modification of the spectrun. Nevertheless, the thermal 
bremsstrahlung function does provide an accurate fit to the spectra 
of most events detected by Mazets ^ energy range 
fron 30 keV to 2 MeV and to the spectrun of the one event detected 
by Gilman et al.® from 2 keV t ■> 2 MeV. Thus, this functional form 
can be considered as representitig the source emission, whatever its 
origin. The deviations fron this spectral form observed at low 
energies in some bursts are thon naturally interpreted as resulting 
fron absorption of this source radiation. 

The observed absorption feature in the X-ray spectrum and its 
variation with time strongly support the model of a neutron star as 
the source of the ganna ray burst. The time variations in the 
total flux on a scale as short as 20 to 40 ms indicate that the 
source size is < 10^ cm. If the burst results from accretion onto 
a neutron star, then the bulk of the radiation probably originates 
at or near the magnetic poles, relatively close to the surface. 
The observations show that the equivalent temperature of the source 
is of the order of 300 keV. .Vs suggested by Mazets al the 
flattening of the spectrum at energies below - 100 keV may be pro- 
duced by absorption at the cyclotron frequency by a cooler plasma 
overlying the source in the presence of a strong magnetic field. 
The fact that the absorption extends over a broad energy range down 
to at least 28 keV indicates that the magnetic field strength within 
the absorbing region must range from <2 x 10^^ gauss to - 10^^ gauss. 
The increase in the energy of the upper edge of the absorption 
feature in intervals 1, 2 and 4 to = 150 keV presumably is a cons- 


equence of higher field strengths in the absorbing region suggesting 
that the source of emission at these tines is closer to the surface 
01 the neutron star. 

The possibility of an emission feature appearing in interval 7 
could indicate that, at this time, the source of emission had moved 
to a greater altitude above the neutron star, similar to the altitude 
of the previous absorption. The absence of any spectral feature in 
the spectrum for interval 8 may result from equal line emission and 
absorption. Alternatively, the emission region may by this time 
have moved to an even greater altitude with much lower magnetic 
fields. Detailed modeling is clearly required before any further 
understanding of these data can be achieved. 

It is interesting to note that two similar though weaker events 
were detected within two days of the April 19 event. One was detected 
two days earlier on 1980 April 17 at 1210 UT with the X-ray detector 
on ISEE-3 and the other two days later on 1980 April 21 at 0308 UT 
with the HXRBS and the GRE on The time profiles of all three 
events show three sharp spikes, and the energy spectra of the April 
19 and 21 events are very similar. Unfortunately, neither event has 
been reported from more than one spacecraft and hence source tri- 
angulations have not been possible. Nevertheless, these three events 
are strongly suggestive of a repetitive nature of the same source. 

In conclusion, if we accept the non-solar origin of this event, 
we have confirmed the results of Mazets et al .^»^ on the existence 
of an absorption feature at X-ray energies in the spectrum of a gamma- 
ray burst. In addition, we have shown variability of the extent of 
the absorption feature on a time scale of 0.5 s during the event and 
detected a possible emission feature at - 100 keV later in the event. 
These results combined with the emission line feature observed by 
Mazets et al .^ at 45 keV in one burst are strongly suggestive of 
the quantization of electron energies in intense magnetic fields in 
the source region. If this interpretation is correct, then the source 
of most gamma-ray bursts must be in the vicinity of neutron stars 
within our galaxy. If this event is, however, of solar or near solar 
origin, then it is a most unusual event and would require reinter- 
pretation on this basis. 

We acknowledge the efforts of the many people who have cont- 
ributed to the success of SMH and the HXRBS. We thank the personnel 
of Computer Sciences Corporation led by Mr. T. 0. Chewning for their 
programming efforts which made the analysis of the data possible. 
We also thank Dr. G. H. Share for many valuable discussions of the 
data. We are grateful to Dr. E. E. Feninore for carrying out Monte 
Carlo simulations of the HXRBS, which were invaluable in determining 
the instrument response as a function of energy, and to Dr. S. R. Kane 
aid Dr. J. M, Ryan for providing X-ray and gamma ray spectra for solar 
events measured on ISEE-3 and with the GRE on SMM. 


10 


REFERENCES 

^•Truoper, J., Pietsh, W. , Reppin, C. , Voges, W., Staubert, R. , and 
Kendziorra, E,, Ap, J. (Letters), 219 , L105 (1978). 

^•Goodnan, N.B., Spac. Sci. Instr. , 1 , 425 (1976). 

^•Mazets, E.P., Golenetskii, W.V., Aptekar, R.L., Gur*yan, Yu. A., and 
Il'inskii, V.N., 1981a, Nature, 378 (1981). 

^•Mazets, E.P., Golenetskii, S.V., Ilinskii, V.N., Gur'yan, Yu. A., 
Aptekar, R.L. , Panov, V.N., Sokolov, I. A., Sokolova, Z.Ya., and 
’Gharitonova, T.V., preprint 719, A. F. IOFFE Physical-Technical Insti- 
tute, Leningrad, submitted to Astrophysics and Sp. Sci. (1981). 

5‘Orwig, L.E., Frost, K.J., and Dennis, B.R., Solar Physics, 25 

(1980). 

^•Fenimore, E.E., Laros, J.G., Klebsadal, R.W., Stockdale, R.E., 
and :>ane, S.R. , preprint, to be submitted to Ap. J. (1981). 

^•Gruber et al . , Ap. J. (Letters), 240 , L127 (1980). 

8‘Gilnan, D., Metzger, A.E., Parker, R.H., Evans, L.G., and Tronbka, 
J.I., Ap. J., 236, 951 (1981). 


^•Matteson, J.L., Ph.D. thesis. University of California at San Diego 
(1972). 

^0‘Share, G.H., Strickman, M.S., Kinzer, R.L., Chupp, E.L., Forrest, 
D.J., Ryan, J.M., Reiger, E., Reppin, C. and Kanbach, G., preprint, 
to be published in the Late Volumes of the 17th International Cosmic 
Ray Conference, Paris, France (1981). 


ON THE THEORY OF GAMMA RAY AMPLIFICATION 
THROUGH STIMtILATED ANNIHILATION RADIATION (GRASAR) 


R. Ramaty, J. M. McKinley*, and F. C. Jones 
Laboratory for High Energy Astrophysics 
NASA/Goddard ^ace Flight Center, Greenbelt, MD 20771 

ABSTRACT 

The theory of photon emission, absorption and scattering in a 
relativistic plasma of positrons, electrons and photons is studied. 
Expressions for the emissivities and absorption coefficients of pair 
annihilation, pair production and Compton scattering are given and 
evaluated numerically. Tbe conditions for negative absorption are 
investigated. In a system of photons and e'*'-e~ pairs, an emission 
line at ~ 0.43 MeV can be produced by grasar action provided that 
the pair chemical potential exceeds 1 MeV. At a ten^erature of ~ 
lO^K this requires a pair density > lO^^cm”^, a value much larger 
than the thermodynamic equilibrivan pair density at this tempera- 
ture. This mechanism could account without a gravitational redshift 
for the observed lines at this energy from gamma ray bursts. 

I. INTRODUCTION 

The emission and absorption of photons in cosmic sources are 
governed by many processes. At teng>eratures of the order 10^ to 
10 ^®K, typical of geunraa ray burst sources, two of the most important 
ones are pair production and annihilation (e^+e~ -»■ Y+Y) and Coit 5 )ton 
and inverse Compton scattering (e+Y J e'+Y')» Arguments based on the 
observed photon intensities of geunma ray bursts and the likely 
distances and sizes of their sources, lead to the conclusion that 
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the source regions of some of the bursts are optically thick. ' 

Photon absorption in yY pair production has been discussed in 
the literature,^ but no calculation has included the effects of the 
stimulation of annihilation or the suppression of pair production 
due to large photon or particle occupation numbers. When these 
stimulation and suppression effects are taken into account, the 
possibility exists for negative absorption.^ The condition for this 
is a population inversion, which in the present context is a pair 
density that exceeds the thermodynamic equilibrium density. 

A recent review of gamma ray burst observations has been given 
by Cline. ^ Of particular interest for the present paper is the 
existence of an emission line seen from several gamma ray bursts in 
the energy range from 0.4 to 0.46 MeV.®“® These lines are probably 
due to e‘^-e“ annihilation radiation. If so, e‘‘'-e“ pairs should be 
present in large numbers in the burst sources, and the sources 
should be sufficiently hot to produce the pairs, but the source 
regions should not be in equilibrium because no lines can then be 

♦Also at Oakland University, Rochester, MI 40063 
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seen. We aim the calculations of the present paper to astrophysical 
sites where such conditions might exist. 

H.KMISSIVITIES AMD ABSORPTIOM COEFPICIEMTS IM A RELATIVISTIC PLASMA 

We consider systems characterized by ten^eratures of the order 
of the electron rest mass energy in which photon-photon collisions 
can produce much larger pair densities than the ambient electron 
densities of tha astrophysical sites of interest. We therefore 
consider only cases in which the electrons and positrons have equal 
densities. As convenient analytical expressions, which allow both 
equilibrium and non-equilibrium situations, we use Bose-Einstein 
distributions for the photons and Fermi-Dirac distributions for the 
pairs. ^ 

We assume equal temperatures for the positrons and electrons, 

® T_» Tj^, but allow the photon temperature, Tyt to differ from 
T^. Since the e*’ and e" densities are equal, n^ » n_ ® n^, these 
particles must also have equal chemical potentials, 

The photon chemical potential. My, is zero for a blackbody 
distribution. We allow non-blackbody photon distributions, but only 
zero or negative values may be assigned to The pair chemical 

potential can be positive, zero or negative. If the pairs 

are in thermodynamic ecpiilibrlum with blackbody photons. 

In terms of these temperatures and chemical potentials, the 
photon and pair densities can be written as^ 

n => (4TrV)"V<a^P n, M) 

where p is the photon or particle momentum and n is the occupation 
number. These are given by 

= {exp[(E^-y^)/kT^l-l}"’ (2) 

for the photons, and by 

n =■ {exp((E -y )/kT ] + 1}"^ (3) 

X XXX 

for the particles, where is the photon energy and E^ is the 
particle total energy (kinetic plus rest mass). 

We proceed to define the photon emissivity and absorption 
coefficient for pair production and annihilation. In particular, we 
are interested in obtaining a correct expression for stimulated 
annihilation which has not been taken into account in previous 
treatments of photon-photon absorption. Stimulated emission has, of 
course, been taken into account for other processes. But we 
cannot use the standard expressions for emissitlvies and absorption 
coefficients because pair production and annihilation do not fit the 
usual pattern in which photons are emitted or absorbed singly and 
the matter has the same form before and after events. We therefore 
proceed as follows: 
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The t.'ansltion.rate in vacuum in either direction between 
photon states in d^p, d^p, and pair states in d\^d^p_ can be 
writt .-n as 

w * d^Pjd^Pgd^p^d^p^ 6 ^(Pj+ 02 -p+-P.)X. (A) 

Here p, and are photon momenta, 6^ and j5_ are momenta of the 
pair, oi. Do, P+ and p_ are the corresponding 4-momenta, and X is 
proportUnaT to the squared matrix element of the interaction, 
summed and averaged over spins and oolarizatlon. 

To take Into account the bath of photons and pairs, we must 
multiply equation (4) by an appropriate expression in the occupation 
numbers incorporating the equilibrium conditions ^ ^ = T and y =u . 
After doing so and integrating over all four momenta, we 'obtain* ^ 
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M Px n.d D_ . 

4ii “11 4ir -h 

1 n,d 0, n^n p, 

4ir fi Att m 


do-i 

DP * 


(5) 


The left-hand side of eguation (*>) is the total pair 
annihilation rate, while the right-hand side is the total pair 
production rate, both influenced by final state occupation 
numbers. The invariant product of the flux factor I and 
differential cross section da/dn is obtained^ from X by integration 
over all final state variables except the angles of one particle: 


[Ic = /p^^dp^/p_^dp_/dp_6^(pj+p2-p^-p_)x. 

The factor in equation (S) is introduced so that each pair of 
photons in either Initial or final state is included lust once. 

To define a photon emissivity and absorption coefficient it is 
necessary to investigate the balance of reactions involving photons 
in a chosen increment d^pj^^ rather than integrating over it. It is 
no longer possible to express the left side in terms o^ the 
annihilation cross section, because the necessary integration of 
eguation (6) has not been performed. But we can interchange the 
order of integration and use equation (7) instead. The required 
balance is then given by 


(6) 

(7) 


(l+nJd^p, dV 

^17?^ ^ aTP^ dn^n+n.(l+n2) 


[Ic 


da-i 
ffn DO 
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■j - 2/dn_^n2[l"n^l (l*n_iric 
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Collecting terms in n^, we obtain 


da-| 

^ PP. 



4 


ori « K A a 


(9) 


4tt (lie) 

where the rate of spontaneous emission, jYy(Ei)» «nd coefficient of 
linear absorption, KyY(Bi^' given by 
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and 
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dp . 

X (El) - /— 3^3/ <sn^.[n2(i-n+)(i-n.) -n^n_( 1+02)1(1^5] • dD 

47T “ft op 

In equation (10) the annihilation emlssivity is expressed in terms 
of the pair production cross section unliVe the approach^^ using the 
annihilation cross section. In the expression for the absorption 
coefficient (equation 11), the first term in the brackets is due to 
absorption by the photon bath while the second term is the 
contribution of induced annihilation. 

Equations (10) and (11) are valid also for nonequilibrium 
situations provided proper nonequilibrium occupation numbers are 
used. In terms of equations (2) and (3), the most general 

nonequilibrium distributions are obtained if T and u fu 

Y t Y * 


With such distributions, the total annihilation and pair production 
rates are not equal and, moreover, Ky can become negative. Ibis 
will happen if an appropriate population inversion takes place. By 
substituting equations (2) and (3) into equation (12) we find that 
for T. » Tv such an inversion occurs if 2u > U • Yn this 


case Kyy is negative for < 2u^-Uy 

For the system to exhibit grasar action, however, it is 
necessary that the total absorption coefficient be negative. For 
the system of photons and pairs that we consider here, the only 
Important process other than pair production and annihilation is 
Cotf^ton scattering. We ignore the weaker processes of 
bremsstrahlung and double Compton scattering. We note, however, 
that synchrotron radiation could potentially be very important, but 
because we are free to choose an arbitrarily low magnetic field 
intensity, we ignore synchrotron absorption in the present paper. 

For Con^ton and inverse Compton scatterings (y,^+e ♦ Y 2 +«' ) we 


proceed in essentially the same way as for pair production and 
annihilation, rising E^ and E 2 for photon energies and p and p' for 
electron and positron momenta, we obtain the Con^ton emlssivity and 
absorption coefficient in the presence of the bath of photons and 
pairs 


cE, 


V=1> 


,d^p 


.da. 


4ir^(flc)^ 


4w^di 
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and 



where the factors of 2 take into account the contributions of both 
electrons and positrons. The Compton emissivity# equation (12), 
represents the scatterings of photons 2 into the element dE^dR^. 

For the absorption coefficient, equation (13), the first term in the 
brackets is due to scattering of photons out of dB^dR^, while the 
second term represents the stimulated scatterings of photons 2 into 
K(E^) can become negative and a necessary condition for 
this is > T^. In our subsequent analysis, however, we shall only 
consider systems with for which Kp is always positive. 

For the numerical evaluations shown mIow we have used the 
expressions of Jauch and Rohrlich^^ for the flux factors and 
differential cross sections. Vte must also express all quantities in 
terms of independent variables of integration. 'Rie detailed 
expressions which we used will be published elsewhere. 

III. NUMERICAL RESULTS 


We have evaluated equations (10) through (13) for various 
choices of and As already indicated, we limit our 

discussion hare to cases with equal pair and photon temperatures, 
■ Ty = T. We allow, however, arbitrary values for and y^. 


We consider first the case of thermodynamic equilibrium, 
Uy “ Ibe emissivities and absorption coefficients for this 


! VT.3ilO*K 
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Fig. 1. Emmissivities and 
absorption coefficients vs. 
photon energy in a system in 
thermodynamic equilibrium at 
3x10^K. The Compton emissivity 
(not shown) la the difference 
between the total emissivltv 
and the annihilation emissivity 
jyy. The Compton absorption 
coefficient (not shown) is the 
difference between the total 
absorption coefficient *nd 
the pair production-absorption 
coefficient Xyy« Tbe photon 
and pair densities in these 
conditions are S.SxIO^^cm**^ and 
2.4x10^^cm“^, respectively. 


case end T >• 3x10^ are shown in *'lgure 1. As can be seen, the 
absorption coefficients are positive at all photon energies. The 
peak of the annihilation emissivity jYY occurs at a higher energy 
than mc^ ~ 0.511 HeV, because the annihilation photons must carry 
away the kinetic energies of the pairs in addition to their rest 0 
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ass energy* 

We turn now to the study of cases with Inverted populations# 

2m. > In this cace# Xyy 1* negative for e < 2m -M , since 
- T^; Kc is positive for all Ey Grasar ait ion 5an’ occur only if 
^t ” Since ty. is proportional to n^^ while the portion 

stimulated annihilation varies as n^^» » sufficiently 
is needed for to exceed K^* l^is implies a 

threshold for m^ which is higher than the threshold required for 
just Kyy to be negative. 

To investigate this threshold we have evaluated - K^y -> K^, 
as a function of m^ for given temperatures and My* We^have c^irried 
out calculations in the teng>erature range 0 < T < 5x10 K# where the 
lower limit corresponds to fully degenerate electrons and 
positrons. We find that the threshold for maser action 1s close to 



and M^ » 1.1 MeV. As can be seen, Kyy is negative for Ey < 2M^ " 


of Kyy due to 
large' density 


2.2 MeV and positive at higher energies. is positive at all 

energies. In the energy range from about 0.25 MeV to 0.7 MeV# “l^yy 



Pig. 2. Bmissivities and 

absorption coefficients vs. 

photon energy in a system at 

3x10^ with blackbody photons 

— and an inverted pair population 

j described by pair chemical 

■j, potential 1.1 MeV. The curves 

■jj have the same significance as 

r in Fig. 1 except that the 

^ Compton absorption coefficient 

S is shown explicitly. 

f Negative values of total and 

“ pair production absorption 

coefficients are represented 

by dashed curves. The photon 

and pair densities in these 

OQ -3 

conditions are 5.5xin'^^cra 
and 7.3x10^^cm"^# respectively. 


(mv) 

exceeds and hence la negative. If the source ie optically 
thick and 10 i* negative over a sufficiently large spatial region# 
then the radiation intensity has a sharp peak at a photon energy at 


which -K id maximal. “Pie value of this eneri 
numerical calculations# is determined prlmarll: 


# ~ 0.43 MeV from the 
from the energy at 


which “Kyy is maximum, shifted somewhat according to the slope of K<;. 
at that energy. Prom Equation (11) %#e can also express Kyy in the 


form 
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Above the graear action threaholdi the first tern* due to stimulated 
emission, Is much larger in magnitude than the second term which Is 
due to absorption. As discussed above, is broadly peaked at 

an energy greater than mc^, reflecting the kinetic energy of the 
annihilating pairs. The division by (from the factor of density 
of states) shifts the peak to an energy somewhat loss than mc^« 

Prom the numerical calculations for other values of T and 
(0<T< SxIO^K, 0.8<y^< 1.2 MeV) we find that the peaks of both 
and -TC^ are in the energy range from about 0,40 to 0.50 MeV. 


IV. ASTROPHYSICAL APPLICATIONS 


The large photon densities expectedl'^ In gamma ray burst 
sources should lead to high pair production and Con^ton opacities. 
Ttie observation®"® of emission lines In the energy range 0.4 to 0.46 
MeV, believed to be due to e'’’-e" annihilation radiation. Is evidence 
that e‘*'-e" pairs do indeed play an lng>ortant role In the physics of 
gamma ray bursts. But It Is not Immediately obvious how a rela- 
tively narrow emission line is produced in a hot and optically thick 
source region. 

Ramaty et al.^®»^® have studied this problem for the transient 
of March 5, 1979^^"^^ from which an emlsslor line was observed® at 
0.43 MeV. 

In their model, the line Is formed In the j optical depth of 
the source region by the annihilation of e^-e" pairs that have been 
cooled by synchrotron radiation prior to their annihilation. The 
observed® upper limit on the width (PWHM <,0.2 MeV), implies^^ a 
temperature less than 3x10®K. it also implies an upper limit on the 
density since even at zero tenq;>erature the line is broadened by the 
motions of the degenerate electrons and positrons. Using eguatlon 
(10) we have evaluated the PWHM of the annihilation emlsslvlty, 

j , as a function of y for T ■ T - 0. 

YY s ± Y 

We find that If FWHM < 0.2MeV, then n^ < 7x10^® cm"^. 

s ~ 

The density can be Independently calculated^® from tlje observed 
line fluence® “ 10 photons cm“^). Let K/{ni) =7.5x10 ' 

cm^sec”^ be the annihilation rate coefficient^^ at 3x10®X, A the 
area of the emitting region. At the time interval in which the 
observed fluence is produced and d “ 5Skpc the distance to the March 
5 source. If the line is formed in a layer of unit optical depth to 
Compton scattering, then 

^ * 2R AAt (4sd^)"\ (15) 

Since R varies as n^^ ** ^ 

With the above numerical values, from Figure 1, and n^ < 
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AAt should oxcasd 1 .5x10^cm^s«c. Ttilg uondltlon is wsll 
satisfisd if the snnihllstlon lins is producsd ovsr ths sntiro 
surfscs of a nautron star# A " lO^^cta^# and during ths ontirs 
ingtulsivs phaas^^ of ths March 5 svsnt# Lt « O.ISssc. Mut ths raodsl 
would facs considsrabls difficultiss if future lasasursiasnts should 
indicats that ths lins is narrowsr than 0.2 MsV# or if AAt should 
turn out, for othsr rsasons, to smallsr than I.SxIO^cm^ssc. 

*^0 advantags of producing ths annihilation lino by grasar 
action is that a narrow lins can Jorm in a hot optically thick 
rsgion. To illustrats this, ««s havs svaluatsd ths photon intsnsity 
psrpsndicular to a slab of thickness L in which ths smissivity and 
absorption cosfficisnt do not dspsnd on pc<^itiont 


I - (j^/K^){1 - sxpt-K^L}}. 


(16) 


Using ths j^'s and K^'s of Figurs 2, %#s show in Figurs 3 ths 
dspsndsncs of I on photon snsrgy ly *hd sl~b thicknsss !>• As can bs 
sssn, grasar action can indssd narrow ths lins in comparison to ths 
thsrraal width'^ of ~ 0.8 MsV in an optically thin Maxws 11-bo It smann 
gas at 3x10^K. 



Fig. 3. Tbs dsvsloprasnt of 
ths annihilation lins with 
inersasing thicknsss of sourcs. 
Tbs systsm is ths same as that 
for Fig. 2. Tbs labels on 
succsssivs maxima indicats ths 
thicknsss involvsd, and ths 
psak snsrgy and FWIM of ths 
lins. 


*:bs psak snsrgy of ths annihilation lins formsd by grasar 
action, is close to the observed peak energies. Tbus, the 
gravitational rsdshift of the line due to the compact object which 
presumably produces tne burst should be quite low, z < 0.1. ibis 
is.pl ies that gamma ray burst sources wi^ observed e'^-s'' emission 
lines cc-uld be objects other than neutron stars, or if they are 
neutron stars, these stars should have msall masoes. 
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Assuming that an invartad layar with paranahara aa in Fiqura 3 
did axiat in tha March 3 burat aourca ragion, tha lina fluanca t Ip 


^ - AAtd“^ /l{l^)dl^ . 

Uaing f.ha raaulta of figura 3 with L “ 10“^ 


(17) 
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'cm and ^ > 10 photona 
cm aguation (17) yialda XAt * 1.2 X 10® cm^aac. By comparing 
with tha minimum hdt daducad abova for radiation producad in tha 
laat optical dapth (AAt > 1.5 x 10^ cm^ aac), wa aaa that not only 
can graaar action produce a narrow lina in a hot region, but that 
tha obaervad lina intanaity and width are conaiatant vrith a much 
amallar aourca and/or a much shorter lina formation time than 
required for tha optically thin case. 


V. COMCI-USIOMS 


Wa have carried out a fully ralvtiviatic treatment of pair 
production and annihilation and Con|>ton and inverse Compton 
scattering in a medium containing photons, positrons, and electrons, 
with equal m* and a” densities. In tha calculation of tha 
amissivitios and absorption coefficients «fa have included tha 
stimulation of transitions caused by tha Bosa-Binstaln nature of tha 
photons and tha suppression of transitions due to electron and 
positron degeneracy. Ws have shown that for systems in 
thermodynamic equilibrium tha calculations lead to an exact balance 
between pair production and pair annihilation. Fbr systems not in 
equilibrium, grasar action is possible. Ws have evaluated, in 
particular, tha absorption coefficient for equal photon and particle 

teng>eratures and positive particle chemical potential 
this example of population inversion, the total absorption 
coefficient can become negative due to the much larger probability 
for stimulated annihilation than for Compton scattering and pair 
production. Graaar action produces a narrow emission line peaked at 
an energy of about 0.43 MeV. ihis energy is lower than the peak of 
the spontaneous annihilation emissivlty, which occurs at energies 
greater than 0.511 MaV. In a bath of blackbody photons (W^ • 0) and 
pairs at the photon temperature, the threshold for grasar 
action is at a 1 MeV corresponding to a pair densities of *- 10^^ 

cm“^ for T » lO^X. A teiq>erature of ~ 5x10^ is needed to produce 
this density in equilibrium with blackbody photona. 

We have applied ouz results to gamma ray bursts, in particular 
to the March 5, 1979 transient from which an emission line at ~ 0.43 
MeV was observed.^ While this line could be produced in a cool skin 
layer, grasar action has the advantage of being capable of 
producing a narrow line from a hot and optically thick source and 
from a source region of relatively small emitting area and short 
duration of line formation. I'Xit if grasar action Is responeible f'-r 
the observed 0.40 to 0.46 Nav emission linej of gamma ray bursts, 
then their sources cannot be a neutron stars of mass larger than 
about 0.6 Mq. the surface of a larger mass neutron star, the 
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gravitational field would shift the line to an energy lower than 
observed. 

There are several difficulties and shortcomings In our 
treatment. We have not shown how the Inversion (v >0) Is 
produced. It could In principle result from the cOolInq of the 
pairs that Is faster than their annihilation, or by a rapid external 
supply of pairs without heating. Cooling by synchrotron emission 
has already been proposed, lb densities that we 

consider here, the required field (B>10l^qauss) seems to lead to l 

synchrotron self-absorption that could quench the grasar action. We 
have Ignored other effects of a strong magnetic field as well, by 
llmitir ’ our calculations to Isotropic distributions and by using ^ 

plane wave functions Instead of Landau functions. This Isotropic 
treatment also does not allow the study of beaming effects which 
should be present In a gamma ray maser. Finally, we have not made 
any attempts to study the spatial and temporal development of a 
system exhibiting grasar action. We expect this development to be 
highly nonlinear. 

As already Indicated, gamma ray burst sources are possible 
astrophysical sites where grasar action could occur. The most 
obvious observational test for this would be the observation of a 
narrow (FWHM « 0.1 MeV) emission line at ~ 0.43 MeV. 
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